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Introduction. The object of this paper is to give an exposition, as ele- 
mentary as possible, of some of those aspects of the theory of the Gamma 
function which are not dealt with in Jensen's "An elementary exposition 
of the theory of the Gamma function. "^ Chapter I presents briefly 

1 Authorized translation, with additional notes, by T. H. Gronwall. These Annals, ser. 2, 
vol. 17 (1916), pp. 124-166. References to paragraphs and footnotes in this paper will be given 
thus: J§2andP. 
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36 T. H. GRONWALL. 

those definitions and theorems on definite integrals which are indispensable 
in the following. In chapter II, the classical applications of the integral 
calculus to the Gamma function are set forth in a form which adheres 
quite closely to the point of view of Jensen's paper, while in chapter III 
the same body of theorems is derived from the definition of r(s) as a 
definite integral.* 

While the paper contains little that is new in subject matter, a con- 
siderable number of proofs have been remodeled, or replaced by new ones. 

CHAPTER I. 

Definite Integrals. 

1. Uniform continuity. Let us consider a function f(xi, Xi, •••, a;„) 
defined for all values of Xi, Xi, • • •■, Xn belonging to a pointset P which 
is bounded (i. e., there exists an A such that | Xi | < A, | a;2 1 < A, • • •, 
I a;„ I < A for all points of P). The function is continuous on the pointset 
P when, for any point X\, Xi, • • -, Xn oi P and any e (arbitrarily small), 
there exists a 6 = 5 (e, Xi, Xi, • • • , Xn) such that for every point xi , Xi', 
• •, x„' of P where 

(1) \ Xi — Xi\ < S, I 3:2' — ^2 1 < 5, • • •, \ Xn' — Xn\ < S 

we have 

(2) 1/(2^1', Xi', ■ •-, Xn) — fixi, Xi, • • •, a;„) I < e. 

The function is uniformly continuous on the pointset P when, for any e, 
a 3 = 6(e) may be chosen independently of xi, Xi, • • •, Xn such that the 
inequality (2) holds for any two points X\, oi>i, • • • , Xn and Xi, x^' , • • • , xj 
of P satisfying (1). The following theorem is of fundamental importance 
for the definition of an integral : 

Theorem I. When f{xi, X2, • • •, Xn) is continuous on a closed^ point- 
set P, this function is also uniformly continuous on P. 

Suppose the theorem to be false ; then there exist an e > 0, an infinite 
sequence of pairs of points on P: Xu, Xi„ • ■ •, Xn^ and XiJ, XiJ , • • •, XnJ 
[y = 1, 2, 3, • • •) such that 

(3) I Xu — a;^ ! < 6„ • • •, i XnJ — a;„J < 6„ 6^ -> as v -» 00, 
and for which 

(4) \f{XxJ, XiJ, ■ •■, Xn,') — f{Xu, Xi„ • • •, Xn,) I S €. 

Let a limiting point of the sequence Xi„ x^,, • ■ ■, Xn, be Xio, Xw, • • ■, Xno', 
this point belongs to P, since P is closed, and /(xi, Xi, • • • , x„) being 

* See the note at the end of this paper. 

2 1, e., a pointset containing all its limiting points. 
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continuous on P, there exists a 5o such that for | Xi — Xio | < 5o, 
I X2 — X20 1 < So, • • •, I x„ — Xno I < 5o we have | /(xi, X2, • • •, x„) — /(xio, 
X20, • • • , x„o) I < «/2. We may now choose v so large that 5^ < 5o/2 and 
I Xu — Xio I < So/2, • • •, I x„^ — x„o I < 5o/2; hence, by (3), | Xi/ — xio| 
S I Xi^ - Xio I 4- I xi/ — xi J < ao/2 + 5^ < ao, • •■•, I x„/ — x„o I < 5o- We 
have therefore | f{xu, X2r, •■•, x„^) — /(xio, X20, • • • , x„o | < e/2, | /(xi/, 
X2/, • • -, x„/) -/(xio, X20, • • •, x„o) I < e/2, whence by addition |/(xj/, 
X2/, • • •, x„/) — fixu, X2y, • • •, x„^) I < e, which contradicts (4), and this 
contradiction proves our theorem. 

2. Upper and lower integrals of a function. We define the interval (a, h), 
where a < 6, as the set of all points x such that a ^ x ^ &. This pointset 
is evidently closed; a and h are called the end points, and the points x 
such that a < X <h the interior points of the interval. The length of the 
interval is & — a. 

A set of values of a real variable | is said to be bounded helow, when 
there exists a finite A such that | > A f or every ^ in the set, and bounded 
above, when ^ < B for every ^. It is a well-known theorem that every 
set of ^'s which is bounded below has a lower hound m such that ^ S to 
for every | but that for any positive e there exists at least one ? in the set 
such that § < vn + e; similarly every set which is bounded above has an 
upper hound M such that ^ ^ ilf for every ^ but ^ > ilf — e for at least 
one § in the set. It is evident that the lower and upper bounds of — | 
are — M and — m respectively. 

In particular, let ^ be the set of values assumed by a function /(x) 
when X varies in the interval (a, h) ; when this set of |'s is bounded both 
above and below, /(x) is said to be hounded in the interval (a, h), and to 
and M are called the lower and upper bounds of /(x) in this interval. 

Now consider a function /(x), bounded in the interval (a, h), and let 
us divide this interval into a number of subintervals (x„_i, x^) (of length 
Xy — Xy-i) by means of the points of division Xo, xi, •••, x„, where 
a = Xo < Xi < X2 < • • • < x„-i < x„ = 6. Denote by M^ the upper, 
and by m^ the lower bound of /(x) in the interval (x^_i, x^), and write 

S = Mi(xi - Xo) + M'2(x2 - xi) + • • • + M„{xr, - x„_i) = 2M„(x, - x,-i), 

S = TOi(Xi — Xo) -t- TO2(X2 — Xi) + • • • + TO„(X„ — X„_i) = Sto^(x^ — X„_i). 

Evidently M S My ^ m^ ^ m, and hence 

(5) M(b -a) ^S ^s ^m{h - a). 

If we consider all possible values of S and s corresponding to all possible 
modes of subdivision of (a, h), S and s are bounded above and below by 
(5) ; hence there exist an upper bound j of s and a lower bound / of *S. 
We shall now prove 
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Theorem II. When the number of subintervals is increased indefinitely 
in such a manner that the length of each tends toward zero, the sums S and s 
tend toward the limits J and j respectively. 

Let us change the sum s into s' by introducing a new division point x' 
between x^_i and x^. Denoting by mj and mj' the lower bounds of /(x) 
in {Xy-i, x') and {x', x^), we obviously have 

s' — s = mjix! — x^_i) + mj'{x, — x') — mX^y — ^v-\) 

= (mj — m^{x' — x^_i) + {mj' — m^{Xy — x'), 

and since ^ m,/ — m^ ^ M — m, OS mj' — m^ ^ M — m, we find 
^ s' — s ^ (M — m){Xy — x^_i) < {M — ih)5 if x^ — x^-i < 5. More 
generally, if s' is obtained from s by the introduction of k new division 
points, we find by introducing them one at a time and successive appli- 
cation of the above inequality that 

(6) S s' - s < k{M - m)b, 

if every subinterval corresponding to s is less than 8. Since — J and 

— S stand in the same relation to — /(x) as j and s to /(x), it is sufficient 

to show that for any e > there exists a 5 such that j ^ s > j — e ior 

any s in which the length of each subinterval is less than 5. Since j is 

the upper bound of all s, it is clear that j ^ s and that there exists an So 

such that So > j — «/2. Let the number of subdivisions corresponding to 

So be k, and consider any s in which the length of each subdivision is less 

than 

. ^ 

^ " 2k{M - m) • 

Finally let s' correspond to the subdivision having as division points all 
those occurring in s and in so] then, by (6), 

s > s' - k(M - m)8 = s' - e/2, 

and applying (6) to So and s', we find s' S So, so that 

S > So - 6/2 > (j - e/2) - €/2 =j - e, 

which proves our theorem. 

The numbers J and j are called the upper and lower integrals of /(x) 
between the limits a and b, and are denoted by 

J = I f{x)dx, j = I f{x)dx. 

Since S S s, it follows from theorem II that J ^ j ot 

(7) fj^'^^'^^ = fj^^)'^^- 
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Let a < c < b and consider the sums S and s for which c is a division 
point, it follows immediately from theorem II that 

/»6 /*c /»6 

I f{x)dx = I f(x)dx + I f(x)dx, 

%Ja %/a \)c 

(8) 

/»o pc n>b 

I f{x)dx = I f{x)dx + I f(x)dxj 

and from (5) and (7) it is seen that 

(9)- m(6 - a) S Cf{x)dx ^ f f{x)dx S MQ) - a). 

Furthermore, when A; is a constant, the sums S and s for lif{x) equal k 
times those for f{x), whence 



(10) 



J/^b /%b /%b /»6 

kf{x)dx = k I f{x)dx, I kf(x)dx = A; I f(x)dx for A; ^ 0, 

a */a «/o »/a 

X6 /»6 /»6 /'ft 

/(;/(x)da; = A; 1 /(a;)dx, ( kf{x)dx = kj f(,x)dx for A; g 0. 



Finally, let/i(x), /2(x), • • -yfnix) all be bounded in (a, 6), and let ilf^ and 
m,, Afi^ and toi^, • • • , Mn^ and m„^ be the upper and lower bounds of 
(fi(x) +/2(x) + ••• +/„(x)), /i(x), •••, fn{x) in (x,_i, X,). Evidently 
M, ^ Jkfu + M2y + • • • + Mny and m^ S mi^ + toj^ + • • • + m„^, whence 

Silf ,(x, — X,_i) S Silfi/x^ — X^i) + • • • 4- ^Mny(x^ — X^_i) 
and 

Sm^Cx^ — Xy-i) ^ Smi^x, — x^i) + • • • + Sm^Xa^^ — x„_i). 

Passing to the Umit, we obtain 

£ ifiix) +Mx) + ■■■ +Ux))dx 

rpb ph 

U{x)dx + I /2(x)dx + • • • + I /„(x)dx, 

J^ (/i(a;) 4-/2(3;) + ••• +/„(x))dx 

fi(x)dx + I fi{x)dx + ■ • • 4- I fn(x)dx. 

3. The definite integral. When the function /(x) is bounded in the 
interval (a, b) and the upper and lower integrals of /(x) between the 
limits a and 6 are equal, the function is said to be integrable in (o, &), and 
the common value of the upper and lower integrals is called the integral 
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of f(x) between the limits a and b and is denoted by 

h 

}{x)dx. 



£ 



The necessary and sufficient condition for integrabiUty is obviously that 
the difference (which is S by (5)) 

S — s = S(M'^ — m^{Xy — Xy-i) 

shall tend toward zero for any one particular mode of subdivision of 
{a, h) when the number of subdivisions is increased indefinitely in such a 
manner that the length of each approaches zero. 

We shall now prove 

Theorem III. A hounded function f(x) is integrable in the interval 
(a, b) when it is possible to find a finite number of subintervals, the total 
length of which is as small as we please, and such that every point of dis- 
continuity of f{x) is interior to one of these subintervals. 

We may note that the hypothesis of this theorem is fulfilled when 
f(x) has a finite number of points of discontinuity in (a, b) . 

To prove our theorem, it suffices, by the remark above, to exhibit a 
subdivision such that the corresponding sums S and s differ by less than 
any assigned e. By hypothesis, we can find a finite number of intervals 

A', of total length less than ktih^-^ — T' ^^^^ t^** every point of dis- 
continuity of f{x) is interior to one of them. Denote by A" the finite 
number of intervals obtained by removing from {a, b) the interior points 
of all intervals A', and consider a subdivision of (a, b) in which the end 
points of all intervals A' are among the division points. Then we may 
write 

S - s = 'E,'{M, - m,){x, - x,_i) 4- ^"{My - m,)(x, - x,_i), 

where 2' and S" refer to the intervals belonging to A' and A" respectively. 
For the former, we have, since M^ — m, ^ M — m, 

S'(M^ - m,){xy - x^i) 

^{M - m)S'(a;, - x^i) < {M - m) ■ 2 (jf'_ ^) = |- 

Regarding S", we observe that A" is a closed pointset on which f{x) is 
continuous, and by theorem I we may assign a 5 so small that 

i/(x')-/(x")i<2(^ 
for any x', x" both belonging to A" and such that \x' — x" \ < 8. Now 
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make every subdivision in A" so small that x^ — x^-i < 8, and let x', x" 
be values of x in the interval (x^_i, x^) for which the continuous function 
/(x) takes its maximum and minimum values M, and m^ respectively. 
Then 

for every subinterval belonging to A", and the total length of these sub- 
intervals not exceeding 6 — a, we have 



Consequently 



^"<2(6-^-(^ -«)=!• 



0^>S — s<24-2=^' 



which proves our proposition. 

Theorem IV. When f{x) is bounded and integrable in (a, b), the same 
is true of | /(x) | . 

It is evident that | f(x) | is bounded when /(x) is bounded. Let M^ 
and m^, MJ and to/ be the upper and lower bounds of /(x) and | /(x) | 
in (x^_i, x^). When My ^ m^ ^ 0, we have MJ = iW„, to/ = to„; when 
S Jlf . ^ My, then Jlf / = — TO,, to/ = — My, and when ikf , g S to,, 
If,' equals the greater of My and — to,, while to,' = 0, so that in all three 
cases ^ MJ — to,' S ikf, — to,. Since S(ikf, — to,)(x, — x,_i) tends 
toward zero, the same is true of '2(My' — to,')(x, — x,_i), so that the 
upper and lower integrals of | /(x) | are equal. 

From (9) it follows that when f{x) is integrable in (a, b) 

(12) m(b - a) ^ f f{x)dx S M{b - a). 

More generally, when three bounded and integrable functions /(x), 
to(x) and M{x) satisfy the inequalities to(x) =/(x) S ilf (x) in (a, 6), 
we have 

m{x)dx ^ I j{x)dx ^ I M{x)dx. 

In fact, since to(x) ^/(x) ^ ilf(x) in (x,_i, x,), the lower bound of to(x) 
in this interval does not exceed that of /(x), and the upper bound of /(x) 
does not exceed that of ili'(x) ; hence the sum s relative to to(x) does not 
exceed that relative to/(x), and the S relative to/(x) does not exceed that 
relative to M{x). Passing to the limit, we obtain (13). 

In particular, theorem IV allows us to take to(x) = — | /(x) | , 
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M{x) = I fix) I , and (13) gives 

(14) I fmdxl^ r\f{x)\dx. 

When a < c < b, (8) becomes for an integrable /(x) 

I f(x)dx = I f(x)dx + I /(x)dx, 
and the repeated application of this formula gives 

fmdx = rf(x)dx + rf(x)dx 

(15) -^^ *^" -^'^ ^.^, ^. 

+ • • • + / /(x)dx + I /(x)da; 

for o < Xi < X2 < • • • < x„_i < 6. 

It is sometimes convenient to consider the case where the upper 
limit in an integral is smaller than the lower one; in this case we introduce 
the definition 

f(x)dx = - I f{x)dx, a <b. 

From (16) it follows that (15) is valid when a, xi, xi, • • •, x„_i, h are in any 
order of magnitude, provided /(x) is integrable in an interval containing 
them all. When h < a, the inequalities (12), (13), (14) must however 
be reversed. 

The inequalities (10) and (11) give when /(x), /i(x), •••, /„(x) are 
integrable 

(17) f kf{x)dx = k i f{x)dx, 

f ifiix) +/2(x) + ••• +/„(x))dx 

= ( /i(x)dx + / /2(x)dx + • • • + I /n(x)dx. 

Theorem V. When /i(x) and f^ix) are hounded and integrable in 
(a, b), the same is true of their -product f{x) = fi{x)fi{x). 

Let M and to, M' and to', M" and to" be the upper and lower bounds 
of /(x),/i(x) and f^ix) in (o, b); furthermore, let My and to„, MJ and m/, 
ilf /' and TO^" be the corresponding bounds in (x^i, Xy). Evidently | /(x) | 
does not exceed the greatest of the quantities M'M", m'm", — M'm", 
— m'M", so that/(x) is bounded in (a, b). To prove the integrabiUty, 
let us first suppose that/i(x) and/2(x) (and consequently /(x)) are never 
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negative in (a, b). Then obviously My ^ MJMy", to„ ^ mjmj' , and 

My — my -^ My'My" " VfiJ ITlJ ' 

= MJ\MJ - my') + my'iMy" - to/') ^ M"{My' - to/) 

+ M'(My" -my"); 
hence 

2(M^ — m^){Xy — Xy-l) S M"'2{My' — my')(Xy — Xy-l) 

+ M'i:{My" - my"){Xy " Xy-l). 

fi(x) and f2{x) being integrable, both sums to the right may be made as 
small as we please by taking the subdivisions sufficiently small, and 
consequently the same is true of the sum to the left, i. e.,f{x) is integrable. 
Second, if /i(x) or/2(x) assumes negative values in (a, b), the product of 
the non-negative functions /i(x) — to' and f2(x) — to" is integrable, and 
by (17) and (18), the same is the case for 

Mx)Mx) = (fi{x) - m')(JM - TO") + to"/i(x) + m'Mx) - m'm". 

Theorem VI. Whenf{x) is integrable in (a, b), and for any x between a 
and b we write 

F{x) = f f{t)dt, x4=a; F{a) = 0, 

then F(x) is continuous in (a, b), and at any point where f{x) is continuous, 
we have 

By (15) we have, when both x and x + h are in (a, b), 

F(x + h) - Fix) = r'mdt - fmdt = r'mdt, 

and if M and to are the upper and lower bounds of fit) in (x, x + h), 
(12) gives 

mh S Fix + h) - Fix) ^ Mh 

(if h is negative, the inequalities must be reversed by (16)), which proves 
the continuity of Fix). When fix) is continuous at x, to and M tend 
toward the same limit fix) as /i -» 0, and consequently 

/(x)=lim^:(^-^M = d_5(5). 
Theorem VII. When f{x) satisfies the conditions of theorem Illy and 
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there exists a function Fi{x) continuous in (a, b) and the derivative of which 
equals f(x) at every point where f{x) is continuous, then 

Fi(x) - Fi(a) = ffmt 

for any x in {a, b). 

Let F{x) be the function defined in theorem VI; then Fi{x) and Fix) 
have the same derivative f{x) in any subinterval of (a, b) in which f(x) is 
continuous, and by the theorem of the mean in the differential calculus, 
Fi{x) = F{x) + const, in the subinterval considered. Now let (xi, Xi) 
and (xs, Xi) be two subintervals of continuity for/(x), separated by the 
subinterval (x2, Xa) in the interior of which f{x) becomes discontinuous. 
Then we have Fi(x) = Fix) + C in (xi, x^) and Fiix) = F{x) + C in 
(xz, Xi), and consequently 

C - C = (Fiixz) - F^(X2)) - (F(x,) - F(X2)). 

By the conditions of theorem III, Xs — X2 may be made as small as we 
please; since Fi(x) is continuous in (a, b) by hypothesis and F{x) by 
theorem VI, the expression to the right may be made as small as we 
please by taking Xs — X2 sufficiently small, and C" — C, being inde- 
pendent of X2 and X3, therefore equals zero. Consequently 

Fiix) = F(x) + C 

throughout the interval (a, b), and since F{a) = 0, we have C = Fi{a) 
and the theorem is proved. 

Theorem VIII. When f{x) , fi(x) andfi{x) are bounded and integrdble 
in (a, b), then 

(19) r f{x)dx = lim S/(?,) (x. - x,_i), 

(20) f fi{x)f2{x)dx = lim S/i(^/)/2(^/')(a;. - x._i), 

where ^v, kJ oii^d ?/' are any points in (x,._i, x^), and x„ — x^_i < 5 /or aiZ v. 
When ^„ is a point of discontinuity of f(x), we may take as the value of f(^J) 
any number between m and M. 

In (19), the sum to the right lies between s and S, and consequently 

tends toward the same limit 1 f{x)dx as these. In (20), we have 

S/l(^/)/2(?/')(^. - X^l) 

= S/i(?/)/2(?/)(x. - x^i) + S/i(^/)(/2(?/') -/2(y))(x. - x.-i); 
the first sum to the right tends toward I /i(x)/2(x)dx as 5 -^ by theorem 
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V and (19), while the absolute value of the second sum does not exceed 

2 |/i(^/) I • |/2(?/') -M^.') I (x. - x^i) 

S 2 |/x(^/) I (MJ' - mj'){x, - x,-{) <{\M'\ + \m'\ )S(ilf/' 

— my"){X, — Xy-l) 

and/2(x) being integrable, the last sum approaches zero as 5 -*0. 

Theorem IX. Integration by parts. Let u(x) and v{x) he continuous 
functions of x such that hounded functions u'{x) and v'{x) exist in {a, b) and 
are continuous and equal to the derivatives of u{x) and v(x) except possibly 
on a pointset which may he enclosed in the interior of a finite number of 
suhintervals of arbitrarily small total length; then 

(21) I u{x)v'{x)dx = u(h)v(b) — u(a)v(a) — I v{x)u'{x)dx. 

Except on the pointset specified, we have 

d{uv) 



uv' -\- vu' = 



dx 



and the integral / {uv' + vu')dx, which exists by theorem V and (18), 

equals u{b)v(b) — u{a)v(a) by theorem VII, whence (21) follows by 
applying (18). 

Theorem X. Integration by substitution. Let f(x) be hounded and 
integrable in {a, b), and (p{t) a function such that <p{a) = a, (pifi) = b, and 
(f>'(t) is continuous and different from zero in (a, /3); then 

(22) rf{x)dx = r f{<p{f)W{t)dt. 

Supposing that (p'{t) > in (a, 0), there exists a A; > such tnat (p'{t) S k 
in (a, 0) . Subdivide (a, /3) at the points a = io, ^i, • • • , <n-i, <n = /3, and 
write Xy = <p{ty). Let M, and m„ be the upper and lower bounds of 
f{(p{t)) in {ty-\, ty) or, which is the same, the upper and lower bounds of 
fix) in {Xy-\, Xy) . Then, by the mean value theorem, 

Xy — Xy-l = <p'(ty'){ty " ty-l) > k{ty — ty-l) > 0, 

where t' is some value in (ty-i, ty), so that 

<Z{My - my)(ty " ty-l) < ^ S (M , " TO,) (X^ " Xy-l) , 

and f{x) being integrable in (a, h) , the sum to the right approaches zero 
when all Xy — Xy-i tend toward zero, which is evidently the case when all 
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ty — ty-i tend toward zero, and the above inequality shows that f{<p{t)) 
is integrable in (a, ;8). The same is true when (p'{t) -^ — k in (a, /3), 
the above inequaUty being then reversed. Since (p'{t) is integrable, the 
product f{<p(t))(p'{t) is integrable by theorem V, and for tj' in (<„_i, U), 
xj = (pity"), we have 

S/(a^.')(x. - x._i) = S/(^(C'))^'(C)(<. - <^i). 

Passing to the limit and using theorem VIII, our proposition follows. 

4. Double and repeated integrals. Let f{x, y) be bounded in the rec- 
tangle a^x^fe, c -^y -^ d, and subdivide this rectangle into the mn 
rectangles x^_i ^ x -^x^, yy-i ^ y ^ yy, where 

a = Xq < Xi < • • • < Xm-i < x„, = b 
and 

C = yo < yi < • • • < 2/n-l < 2/n = d. 

Denote by M and m the upper and lower bounds of /(x, y) in the rectangle 
a^x^b, c^y^d, and by M ^^ and m^^ the corresponding bounds in 
the rectangle x„_i ^ x ^ x^, yy-i ^ y ^ yy, and consider the two sums 

m n 

S = Y, T.M^y(x^ - x^-i)(yy - yy-i), 

m n 

s = S Uniy.yix^ - x^-i)(yy - yy-i). 

(i-1 y=l 

Since evidently M S M^^ S to^„ S to, we have 

(23) M(& - a)(d - c) m S ^s ^ m{b - a){d - c), 

so that for all possible modes of subdivision of the rectangle a ^ x ^ fe, 
c ^ y ^ d, S and s are bounded above and below; hence there exist an 
upper bound j of s and a lower bound J of S. In exactly the same way 
as we proved theorem II, we find 

Theorem XI. When the number of subrectangles is increased indefinitely 
in such a manner that the sides x^ — x^-i and yy — j/^-i of each tend toward 
zero, the sums S and s tend toward the limits J and j respectively. The 
numbers J andj are called the upper and lower integrals of /(x, y) over the 
rectangle a^x^b, c^y^d, and are denoted by 



J = J J /(^' y)dxdy, j = j J fix, y)dxdy. 

When/(x, y) is bounded in the rectangle, and the upper and lower integrals 
of /(x, y) over the rectangle are equal, the function is said to be integrable 
in the rectangle, and the common value of J and j is called the (double) 



THE. THEORY OF THE GAMMA FUNCTION. 47 

integral of f{x, y) and denoted by 

The necessary and sufficient condition for integrability is obviously that 
the difference (which is S by (23)) 

TO n 

S - s = Y.Yj {M^, - m^;){x^ - x^-i){y, - j/^i) 

shall tend toward zero for some particular mode of subdivision of the 
rectangle when the number of subrectangles is increased indefinitely in 
such a manner that the lengths of all their sides approach zero. 

The following theorem is proved in exactly the same way as theorem 
III: 

Theorem XII. A hounded junction f(x, y) is integrable in the rectangle 
a^x^b, c^y^d, when it is possible to find a finite number of sub- 
rectangles, the total area of which is as small as we please, and such that 
every point of discontinuity of fix, y) is interior to one of these subrectangles. 

For any value of x in (a, b), f{x, y) is bounded when c S y ^ d, and 

Xd pd 

f(x, y)dy and I f{x, y)dy 

exist. Denote by _ 

Fix) = £f{x, y)dy 

any number such that 

fix, y)dy S j fix, y)dy S j fix, y)dy, 

r*d 

so that, in particular, Fix) = I fix, y)dy for every value of x for which 
fix, y) is integrable in respect to y between the limits c and d. Then the 
integral (which we suppose for the moment to exist) I Fix)dx is called 
the repeated integral of fix, y) in respect to y and x and is denoted by 

Xb /%d 

dxj fix, y)dy. 

Similarly, we may define the repeated integral in respect to x and y, 

dy I fix, y)dx = I Giy)dy, 

where Giy) = I fix, y)dx. We then have 
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Theorem XIII. When f(x, y) is hounded and integrdble in the rectangle 
a^x^b, c^y^d, then both repeated integrals of f{x, y) exist and are 
equal to the double integral: 

dx I fix, y)dy = \ dy \ f{x, y)dx = I ( f{x, y)dxdy. 

Since m^, ^f(x, y) ^ M^^ for x^-i ^ x ^ x^, y,-i ^ y ^ y^, it follows 
from (9) that 

fix, y)dy ^ I fix, y)dy ^ M^Xv. ~ y.-d; 
adding these inequalities for i- = 1, 2, • • •, n and applying (8), we find 

Hm^.iy, - y^-i) ^ I fix, y)dy ^ | fix, y)dy ^ Y.M,,Xy. - Vv-i) 

for any x in (x^_i, x^), and consequently the lower bound m^ and the 
upper bound ilf^ of Fix) = I fix, y)dy for x^_i S x S x^ satisfy the 
inequalities 

TO n 

Multiplication by x^ — x,^-! and summation in respect to ^ now gives 

m n m 

Z Jlm,,,iXi. - x,,-i)iy, - y,-i) ^ J2 rn,,ix,, - x^_i) 

^ Z -M'^Cx^ - x^_i) 

m n 

S E Z-M"^.(x^. - x^_])(2/, - 2/,-]). 

But as m and n tend toward infinity while all x^ — x^_i and y„ — y^-i 
tend toward zero, the two outside members of this inequaUty approach 

the same limit ( ( /(x, y)dxdy; consequently the two inner members 
must also approach this same limit, that is, I Fix)dx exists and equals 

Ja 

the double integral. The same proof obviously applies also to the other 
repeated integral. 

df(x Vi 
Theorem XIV. When fix, y) and ^ o,re bounded for x in ia, b) 

and y in (c, d) and 
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exists for any y in (c, d), 



rdf{x,t) 



dt 



dt 



exists for any x in (a, b) and any y in (c, d), and finally I f{x, c)dx exists, 

Ja 

then 

for every value of y for which the last integral is a continuous function of y. 
It may be noted that all the conditions of the theorem are satisfied 
when/(x, y) and df{x, y)/dy are continuous for x in (a, b) and y in (c, d). 
For the proof we observe that by theorem XIII 

and by theorem VII 

f'^-^dt=f{x,y)-f{x,c); 

hence, by theorem XIII, I (/(x, y) —fix, c))dx exists, and from the 
existence of I fix, c)dx and (18) we conclude that 

I f{x, y)dx = I fix, c)dx + I (/(x, y) - fix, c))dx. 
Consequently (24) gives 

(fix, y)dx - ffix, c)dx = ("dt r^^dx, 
and differentiating in respect to 2/, theorem VI gives 

for any value of y for which the integral to the right is a continuous 
function of y. 

The following proposition is proved in exactly the same way as 
theorem IV: 

Theorem XV. When fix, y) is bounded and integrable in the rectangle 
a ^ X ^ b, c ^ y S d, the same is true of | /(x, y) \ . 

Formulas (13) and (14) are also immediately extended to double 
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integrals: 

(25) I I mix, y)dxdy ^ / I f(x, y)dxdy ^ I I M(x, y)dxdy 

when the three bounded and integrable functions f{x, y), m{x, y) and 
M{x, y) satisfy the inequaUties m{x, y) ^ f{x, y) S M{x, y) in the rec- 
tangle of integration; andfor TO(a;, y) = — |/(x, y) \ , M{x, y) = |/(x, y) \ , 

(26) I f f /(x, 2/)dx% I ^ r r I /(x, 2/) I dxdy. 
Also (17) and (18) extend immediately to double integrals: 

kf(x, y)dxdy = h \ 1 /(x, y)dxdy, 

%Ja vc 

I (/i(a;, y) + /2(a;, 2/) + • • • + /„(x, y))dxdy 

ph pd /»6 pd 

(28) = /i(x, y)dxdy + /2(x, y)dxd2/ + • • • 

X6 /^d 
J fn(x,y)dxdy. 

In (25)-(28), any double integral may evidently be replaced by either 
of the repeated integrals. Finally, theorem VII extends immediately to 
double and repeated integrals. 

5. Infinite integrals. In the preceding paragraphs, we have assumed a 
and b finite and/(x) bounded in (a, b). When one of these conditions is 
not fulfilled, we may define, by a limiting process, what is called an infinite 
integral (the term is chosen in analogy to "infinite series"). 

First assume /(x) to be bounded and integrable in (a, b) for any b > a. 
The integral between the limits a and <x> is then defined as 

I f(x)dx = lim I f{x)dx, 

whenever this limit exists. Similarly we may define 

}{x)dx = lim I f{x)dx, I f(x)dx = hm I f(x)dx. 

»->oo 

Next let us assume that /(x) is bounded and integrable in (a, b — e) 
for an arbitrarily small positive e, but tends toward + °° or — oo (or 
oscillates between these bounds) as x -» ft. The integral between the 
limits a and b is then defined as 



I f{x)dx = lim I /(x)da 
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whenever this limit exists. Further, when f(x) ceases to be bounded as 
a; -> o, we define 

f{x)dx = lim I f(x)dx, 

and when c is a point between a and b in the neighborhood of which f(x) 
is not bounded 

I f{x)dx = lim I f(x)dx + lim I f(x)dx. 

Similarly we may define infinite repeated integrals; for instance, when 
f(x, y) is bounded and integrable for a + « S x ^ 6, c + e' ^ y ^ d, 
where e and e' are as small, 6 and d as large as we please, we have the 
definitions 

I dx ( /(a;, y)dy = lim I do; lim I /(a;, 2/)d2/ , 
I dy I /(x, 2/)da; = hm I % Um I /(x, 2/)dx . 






In the following, we shall consider only integrals of the type 
f(x)dx where f(x) may cease to be bounded at x = a, and the corre- 



/ 

I/O 

spending case for repeated integrals. The theorems we shall prove hold, 
however, in the other cases; the necessary modifications of the proofs are 
obvious. 

A convenient means of establishing the existence of an infinite integral 
consists in comparing it to another with positive integrand. 

Theorem XVI. Letf(x) be bounded and integrable when a + e ^ x ^b, 
where e is as small and b as large as we please, and let | /(x) | ^ M(x) in this 

interval. If ( M{x)dx exists, and is bounded when e -» and 6 -> oo , 
f{x)dx exists, and 

/(x)dx g I M{x)dx. 

Since M(x) S 0, | M(x)dx increases (or remains constant) when e 
decreases or 6 increases, and being bounded, this integral must therefore 

XOO 
M{x)dx exists. 
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Consequently 



lim f M{x)dx= 0, lim f M(x)dx = 0, 

e'-»0 &'-»oo 

and by (14) and (12) 

I r^' f{x)dx\ ^ P' Mix)dx, f fix)dx\ ^ f Mix)dx, 

BO that 

f{x)dx = 0, lim I f{x)dx = 0, 

Xoo 
/(a;) da;. 

Corollary. TF/ien |/(x) | < K{x — a)"", where a <l,for a < x < Xi, 

/(x)dx exists. 

a 

In fact, when e and e' are so small that a + e < a + e' < a;i, we have 
I f*' f{x)dx\^ f' Kix - a)-'dx = 3-§-(A" - e'-°) -*0, 
and for 6' > 6 > X2, 

if /(x)dx| sf X.-^dx = ^^(^-^-j,,) -.0. 

When I f{x)dx exists, then, by the definition of an infinite integral, 

Jr»ii /»oo 

f(x)dx and I f{x)dx also exist, and (15) becomes 

/•oo /»J:i /»*2 /'•'n-l /^"O 

(29) fix)dx = f{x)dx + I /(x)dx + • • • + /(x)dx + j{x)dx 

Ja Ja 'Jll «^^n-l 'J^n-l 

for a < xi < xs < • • • < x„_i. The equation of definition (16) becomes 

(30) r /(x)dx = - r /(x)dx, 

and (17) and (18) also extend immediately to the present case: 

A;/(x)dx = k I /(x)dx, 
/" (/i(x) +/2(x) + • • • + /„(x))dx 

Mx)dx + Mx)dx + • • • + /„(x)dx. 
By a passage to the Umit in theorem VII, we find 
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Theorem XVII. When f(x) satisfies the conditions of theorem III in 
(a + e, h) where e is as small and b as large as we please, and there exists 
a function Fi(x) continuous in (a + e, b) and the derivative of which equals 
fix) at every point where f(x) is continuous, when finally 

lim Fi(a + e) = Fi(a) 

Xoo 
f{x)dx exists and 

Fi(«) -Fi(a) = rf{x)dx. 

Similarly, it is seen that when the conditions of theorem IX are fulfilled 
in (o + e, 6), where e is as small and 6 as large as we please, and three 

u{x)v' {x)dx, I v{x)u'{x)dx, u{a)v{a) = 

*) a 

hm u{a + i)v{s + e) and «(oo)t;(oo) = lim «(&)«;(&) exist, then the fourth 

e-> S-»go 

exists and 

(33) I u{x)v'{x)dx = ■w(oo)j;(oo) — uia)v{a) — | v{x)u'{x)dx. 

Finally, it is readily seen how theorem X, on integration by substitution, 
may be extended to the present case. 

A convergent infinite series Mo + 'Wi + • • • + «„ + • • • may be written 
as an infinite integral by making /(x) = «„ f or n ^ x < n + 1 . In the 
interval (0, 6), where n + l^&<n + 2,/(a;) has the only discontinuities 
X = 1, 2, • • -, n + 1 and is integrable by theorem III, and (15) gives 

I f{x)dx = Z) I fix)dx + \ f(x)dx = Z^M^ + (& — n — !)«„+!. 

The series being convergent, Un+i -> as n ^ <» , so that 
\ (b — n — l)u„+i I < I Un+i I -^ 0. 

In the following, the concept of uniform approach to a limit (or uniform 
convergence) is of fundamental importance. Let f(x, 8, w) depend on the 
parameters 5 and w; then/(a;, 5, w) tends toward the limit /(x) uniformly in 
(a, 6) as 8 ^ and w -> oo when, for any e > as small as we please, 
there exist a 6o = 8o{e) > and an wo = wo(€) > independent of x such 
that for every x in (a, 6) 

I /(x, 5, w) - /(x) I < € for < 5 < 8o, CO > coo. 
In particular, let f{x, 8, co) be independent of 6 and make 

f(x, w) = Z^Xa;) 
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for n ^ w < n + 1; we then obtain the famiUar definition of uniform 
convergence of a series. 

Theokbm XVIII. When for a ^ x ^b, fix, 6, w) is hounded and 
integrahle in respect to x and tends uniformly toward f{x) as 8 -^ and 

w -> 00, then I f{x)dx exists and 

rb /*h 

f{x, 6, o})dx = I f{x)dx. 






For any e > 0, there exist, by hypothesis, a 5o = SoCe) and an coo = wo(€) 
such that f{x, 5, w) — e ^ f{x) ^ f{x, 5, co) + e in (a, 5) for < 5 < 5o, 
w > Wo- Hence f{x) is bounded in (a, h) so that its upper and lower 
integrals exist; the application of (11) and (9) to the above inequality 
gives _ _ 

Xb r*b 

f(x)dx ^ I fix, 5, <S)dx + €(5 - a), 
*J a 

I fix)dx ^ I /(x, 6, w)dx — €(& — a). 

tJa ^ a 

Since /(x, h, (a)_ is integrable, its upper and lower integrals are equal, 
fix)dx — I fix)dx S 2€(6 - a); but fix) being inde- 

pendent of «, it follows that | fix)dx exists, and the preceding inequalities 
become 

fix)dx - I /(a;, 5, co)dx ^ «(& - a) 

for < 6 < 5o, w > coo, which proves our proposition. 

fix, y)dx is said to exist uniformly for y in (c, d) when 
/(a;, 2/) is bounded and integrable in respect to x in (o + e, h) and the 
integral I fix, y)dx tends toward I fix, y)dx uniformly for y in (c, d), 
that is, when 

fix, y)dx = 0, lim fix, y)dx = 

e'-»0 »'-» M 

uniformly for y in (c, d). A convenient test for uniform existence is 
contained in 

Theorem XIX. When fix, y) is hounded and integrahle in respect to x 
in (a + e, h), while y is in (c, d), and when, for x > a and y in (c, d), 
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uni- 



I M{x)dx exists, then 1 f{x, y)dx exists 

formly for y in (c, d). 

In fact, the proof of theorem XVI then shows that both the integrals 

fix, y)dx and I f{x, y)dx 

approach zero uniformly as «, e' -> and 6, &'-><» . 
We furthermore have 
Theobem XX. When f{x, y) is hounded and doubly integrable for 

a + e-^x'^byC'^y^d, and I f(x, y)dx exists uniformly for y in (c, d), 

Ja 

then 

dy I fix, y)dx = \ dx \ fix, y)dy. 

By theorem XIII, we have 

\ dy I fix, y)dx = i dx j fix, y)dy, 

and since I f{Xj y)dx exists uniformly for c ^ 2/ = ^j theorem XVIII 

/'*d /^oo f*d /»6 

shows that | dy 1 /(x, y)dx exists and equals lim I dy I /(x, ?/)dx. 

/•ft /»d 

The preceding equation then shows that lim I dx 1 fix, y)dy exists 

e-»0 Ja+e Jc 

and equals | dy ( fix, y)dx. But the last limit, by the definition of an 
infinite integral, equals I dx I fix, y)dy, whence our theorem. 

%Ja «/c 

Making /(x, y) = «„(2/) forn ^ x < n + 1, where all Uniy) are bounded 
and integrable for c ^ y ^ d, the uniform existence of I fix, y)dy is 

Jo 

00 

equivalent to the uniform convergence of Y^Uniy) in (c, d), whence 



writing x, a, h for y, c, d, we find the 

00 

CoEOLLARY. When TjUnix) converges uniformly in ia, b), every u„ix) 



bein^ bounded and integrable, then 

y^6 00 CO /»6 

I ^Un{x)dx = ^ I Un{x)dx. 
•Ja t/a 
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Theorem XXI. When f{x, y) and df{x, y)ldy are continuous for 
a + e^x^b, c^y^d, I f(x, c)dx exists and I {df(x, y)/dy)dx 

exists uniformly for y in (c, d), then I f{x, y)dx exists uniformly for y in 
(c, d), and 

vyi f(^'y)d^ = i -^r~^^- 

All conditions of theorem XIV being fulfilled for x in (o + e, b) and y in 
(c, d), we have 

r f(x, y)dx = f fix, c)dx + {"dt C ^%^dx. 

Since the first integral to the right is independent of y and tends toward 
f{Xj c)dx as € -> 0, & -> 00, and since I dt 1 (c'/(a;, t)ldt)dx exists uni- 

Xco 
{dj{x, y)ldy)dx and theorem 

XVIII, it follows that I /(x, 2/)<^^ exists uniformly for y in (c, d) and 

£f{x, y)dx = £f{x, c)dx + £dt£°°^-^^^dx. 

From the continuity of 6i/(a;, t/)/^?/ for a + e ^ x ^ fe, c ^ y ^ d, and 

the uniform existence of I (3/(x, y)ldy)dx for y in (c, d), it follows at 

once that the latter integral is a continuous function of y in (c, d). Con- 
sequently, differentiating the last equation in respect to y and applying 
theorem VI, we obtain our proposition. 

Making /(x, y) = Un{y) for » ^ x < n + 1, and replacing y, c, d by 
X, a, b, we have the 

00 

Corollary. When ^(dun(x)ldx) converges uniformly in (a, 6), every 



00 CO 

{dun{x)ldx) being continuous, and X)«n(«) converges, then 23«„(x) converges 



uniformly in (a, b), and 

^2:«„(x)=E-^— . 

Regarding infinite repeated integrals, we shall now prove 
Theorem XXII. When /(x, y) is bounded and integrable for a + e ^ x 
^ b, c + e' ^ y '^ d, where e, «' are as small and b, d as large as we please, 
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and ONE of the infinite repeated integrals I dx \ |/(x, y) \ dy and 

I ^2/ I |/(^)2/)|^^ exists, then both the infinite repeated integrals of 
/(x, ?/) exis< and are equal: 

I dx I f{x, y)dy = j dy j fix, y)dx. 

I /(^) y)dxdy we infer, by theorem XV, 
I l/(^) y) I da:^?/ a^d hence by theorem XIII that 

-. + € tA+e' 

I dx I l/(x, 2/) I dy = I dy I |/(x, 2/) | dx. 

When in either of these repeated integrals e and e' are decreased, & and d 
increased, the integral cannot decrease since |/(x, y) \ is never negative. 

dx I |/(x, 2/) I dy, we have 
d^ I |/(^. 2/) Ml/ = I dx \ |/(x, 2/) I dy, whence 

dy I \f{x, y) I dx ^ I dx \ |/(x, y) | dy, 

SO that the integral to the left is bounded for all values of «, e', h and d 
considered. Since this integral does not decrease as e, «' decrease to zero 
and 6, d increase indefinitely, it therefore approaches a limit not exceeding 

% I I /(s^) y) I dx exists and 

c ^ a 

I dy |/(x, 2/) I dx S I dx \ |/(x, y) | dy. 

Having established the existence of I dy \ \ f{x, y) \ dx, we apply to it 

the same argument as was applied to the other infinite repeated integral, 
and thereby establish the inequality 

dx |/(x, y) I dy ^ I dy \ |/(x, y) | dx. 

Combining the two inequalities, we find 

dx I |/(x, y) I dy = I dy I |/(x, y) | dx. 
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Furthermore, the integral I I (|/(x, y) \ — f(x, y))dxdy exists, being 

the difference of the double integrals of |/(x, y) \ and f{x, y), and from 
S I f{x, y) I - J{x, 2/) g 2 I j{x, y) \ we see that 

Q^C dxf (|/(x, y) I - fix, y))dy S 2 T dx f" \f(x, y) \ dy. 

Since the integral to the left cannot decrease as «, «' decrease to zero and 

b, d increase indefinitely, it follows that | dx j (\f(x, y) | — /(x, y))dy 

exists; we may now show exactly as before that the other infinite repeated 
integral exists and that 

r dx r (\f(x, y) I -fix, y))dy = f dy f (|/(x, y) \ -fix, y))dx. 

Our theorem now follows by subtracting this equation from that relative 
to \fix, y) I . 

Replacing fix, y) by Uniy) for » S a; < » + 1 and writing x, a for 
y, c we have the 

Corollary. When every Unix) is finite and integrable in (a + «, b) 

/»0O 00 CO /»C0 

and either of the expressions j 1^\unix)\ dx and J2 i \unix)\ dx exists, 

Ja Ja 

then both the corresponding expressions in Unix) exist, and 

J<*00 00 00 /»00 

Y^Unix)dx = 53 I Unix)dx. 
a J a 

6. Integrals of complex functions of a real variable. In the following two 
chapters, we shall make extensive use of such integrals, which are defined 
in the following manner: let/(x) = giix) + ig^ix), where giix) and g^ix) 
are real functions of x, bounded and integrable for a ^ x S 6. Then, by 
definition, 

fix)dx = I giix)dx + i j g2ix)dx, 

and the double and repeated integrals oi fix, y) = giix, y) + igiix, y), as 
well as the various kinds of infinite integrals, are defined in the same 
manner by decomposition into their real and imaginary parts. 

Theorems II and XI, as well as formulas (7)-(13), apply to real func- 
tions only; but all the remaining theorems and formulas may be extended 
to the present case. In fact, the proof of theorem I subsists unchanged, 
since it makes no use of the reality of /(xi, • • •, x„). Theorem III is 
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immediately extended by decomposition of f(x) into its real and complex 
parts. In theorem IV, let M^, rriy, Mi„, mu, M^^ nhy be the upper and 
lower bounds of ] /(x) ], ] gi(x) \ and ] gi{x) | in {Xy-i, Xy), and x', x" any 
two points in this interval. We have 

\f(x')\-\fix")\^\f(x')-fix")\ 

= 1 9iix') - g,{x") + i{g,{x') - g,{x")) \ 

^ I g^{x') - g^{x") I + \g,{x') - g,(x") \ 

^ (Miy - miy) + {M2y - nhy), 

and since for suitably chosen values of x' and x" \ f{x') \ will differ from 
My, and | J{x") \ from niy, as little as we please, it follows that 

My - rriy -^ {Miy — mu) + {M2y - nhy), 

and consequently the S — s relative to | f{x) \ may be made as small as 
we please, since it is less than the sum of the S — s relative to | gi(x) \ 
and that relative to | g2{x) \ , both of which may be made as small as we 
please by theorem IV. Theorems V-X, XII-XIV are immediately 
extended by decomposition (it being observed that in X (f>(t) is a real 
function of the real variable t), XV is extended in the same way as IV, and 
the proof of XVI remains unchanged, since (14) subsists in the complex 
case. In fact, with the notations of theorem VIII, we have, since 
Xy — Xy-1 > 0, for/(x) complex 

1 S/(L)(x. - x._i) I ^ S 1/(L) 1 (x. - x_i); 

as we have seen that theorem VIII extends to the complex case, it follows 



that for 8 -» 0, the left side approaches 



I f(x)dx 



and the right, 



I 1 /(x) I dx, whence we obtain (14). In exactly the same way, we extend 

(26), while (15), (16), (17), where A; may now be a complex constant, 
(18), (27), (28), (32), (33) are extended immediately by decomposition. 
Theorems XVII and XVIII extend by decomposition, the proof of XIX 
remains unchanged, XX and XXI (with corollaries) are again extended 
by decomposition, while in XXII (and corollary) we first conclude from 
I gxix, y)\^\ /(x, y) \ and | g^ix, y)\-^\ fix, y) \ that the same infinite 
repeated integrals of | gi(x, y) \ and \gi(x, y) \ exist as that of |/(x, y) | 
included in the hypothesis; the extension is then finished by decomposition. 

We shall finally prove the following 

Theoeem XXIII. When, for an e > as small and ab as large as we 
please, the function f{s, t) of two complex variables s and t is holomorphic for 
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a + e ^ t ^b and \ s — So\ S r, where r is independent of e and b, and 

/»oo 

when finally I /(s, t)dt converges uniformly for \ s — So\ < r, then this in- 

tegral defines a function of s which is holomorphic for \s — So\ <r. 

Write s — So = X + iy and /(s, t) = gi{x, y, t) + ig-iix, y, t), and 
denote by P the pointset defined hy a + e ^ t ^b, x^ + y^ ^ r^; since 
f(s, t) is holomorphic on P, gi and grj are continuous on P and this pointset 
being closed, gi and g2 are uniformly continuous on P (theorem I). Since 
^1 and g2 are continuous, we may take the whole interval (a + «, b) as A" 
in theorem III, and on account of the uniform continuity on P, we may 
choose, for any e' > 0, a 5 independent of x and y such that >S and s differ 
by less than e' when all subdivisions of (a + «, b) are less than 5. Conse- 
quently, as 5 -> 0, the sum 

2/(s, T,)(<, - <,_i) = SflfiCx, y, r,){t, - <,_i) + i1,g2{x, y, r,){t, - <,_i), 

formed according to theorem VIII, tends toward I /(s, O^^' uniformly 

for I s — So I = r. But each /(s, t„) is expansible in a power series in 
s — So convergent for | s — So | ^ r, and therefore also S/(s, t^) {t, — i^_i) 
which contains a finite number of such terms. Let 8i, h, • • • be a sequence 
of positive numbers such that 5„ -> as « -» oo , and let 

*„(s; e, b) = S/(s, T,)(i, - <,_i) 

for a certain mode of subdivision where all t^ — ty^i < 5„; since f>„(s; e, 6) 
is a power series in s — So convergent for | s — So | ^ r and since 

$„(s; «, b) tends uniformly toward | /(s, <)d<, we conclude from Weier- 

strass' theorem on sequences of analytic functions that I /(s, 0^^ is 

expansible in a power series in s — So convergent for | s — So | ^ r. 
Now choose sequences ti, e^, • • ■ and &i, 62, • • • such that e„ -> and 

/(s, Oo!' converges uniformly for | s — So | ^ r, 
the expressions I /(s, t)dt, which are power series in s — So convergent 

/(s, t)dt which 

is therefore, by a second application of Weierstrass' theorem, a power 
series convergent for | s — So | < r. 

7. Fourier series. Let/(a;) be bounded and integrable for « ^ a; ^ 1 — e, 

where « is as small as we please, and suppose that I |/(x) | dx exists. 
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Since | /(x) sin 2mrx \ ^ | /(x) | and | f{x) cos 2nvx \ S | f{x) | , it follows 
from theorem XVI that the integrals 

(34) a„ = 2 I f{x) cos 2nTxdx, 6„ = 2 I f(x) sin 2mrxdx 

t/O t/0 

both exist. Then the Fourier series of f{x) is defined as 

(35) "o + Z^ (a» cos 2n7rx + b„ sin 2n7rx). 

In regard to the convergence of this series, we shall need in the following 
only this very special theorem: 

When /(x) is holomorphic for < x < 1 (but not necessarily at the end 

points of the interval) and I |/(x) | dx exists, the Fourier series of f(x) 

Jo 
converges to the sum f{x) uniformly for e ^ x ^ 1 — e, where e > hut as 
small as we please. 

Using the integration variable t instead of x in (34), we find 

o„ cos 2n7rx + bn sin 2n7rx 

= 2 I f(t) (cos 2mrt cos 2nvx + sin 2nirt sin 2mrx) dt 
Jo 



= 2 \ f{t) cos 2mv{t - x)dt, 

Jo 



and since we have' 

1 I o -T^ n r4 \ sin (2n + l)7r(< - x) 

1 + 2 ^ cos 2vK{t — x) = ^^ ., \ 

y=i sm ir{t — x) 

it follows that 



Oo 

2 



+ 23 {a, cos 2j'7rx + by sin 2vvx) 



v=l 



(36) = j[ /(O [ 1 + 2 i: cos 2vi:{t - x) j di 

^ r- sin(2n;Wx) ^^_ 
Jo sm irit — x) 



' Proof; in the finite geometric series 



yto e^»i - I e»i-e-»* 

_ e(o+(2n+l)S)t _ g(a-ff)i 

2i sin /3 
compare the real and imaginary parts on both sides: 

2 cos (« + 2.« = «m(a + (2« + m_-sin(a-ff) 

2 sin (a + 2.0) = '^°«(«+^)-c°s„( <' + (2« + lW 
v=o ^ < 1^' 2 sm /3 

From the first of these, we obtain the formula in the text by making a = 0, /? = 7r(< — x). 



62 T. H. GRONWALL. 



In the special case f(x) = 1, (34) gives ao = 2, a, = b^ = {v > 0), and 
(36) becomes 

'' sin (2n + l)7r(< - x) 
^0 sin ir(i — x) ' 



-r 



multiplying this by/(x) and subtracting from (36), it is seen that 
2 



+ S (dy COS 2;'7ra; + b^ sin 2fxx) — f(x) 



(37) = r (/(O _ /(^)) sin(2n + l)x(t-a:) ^^ 

= /l + •72 + Ji, 

where Ji, Jz and Jz are the integrals corresponding to the decomposition 
of the interval (0, 1) into the three intervals (0, 8), (5, 1 — 5) and (1 — 5, 1). 
Our theorem will then be proved if we show that after assigning an 
arbitrarily small e, a 5 may be chosen sufficiently small and an no suffi- 
ciently large, both independent of x, such that \Ji\ + | J2 1 + | J3 1 will 
be as small as we please for n > no and every x in («, 1 — e). Assume 
< 5 < e/2; for e S x S 1 - « and 0^<S5 or 1-5^<^1, we 
then have e/2 < | < — x | < 1 — (€/2) and consequently | sin Tr{t — x)\ 
> sin (7re/2); since \f(t) — f{x) \ S \f{t) \ + \f{x) |, the expression under 
the integral sign in Ji and J 3 will be less than ( | /(<) | + | f(x) \ )/sin ir(e/2), 
whence 



1 Ji 1 + I /a I < 



'^£(\m\ + \m\)dt 



sin (7re/2) 



+f\i\m\ + \f{x)\)dt^. 



As f(x) is holomorphic in («, 1 — e), there exists a constant M such that 
I f{x) I < ikf for all x in («, 1 — e) , whence 

Since I 1/(0 | dt exists by hypothesis, the last two integrals may be 

Jo 
made as small as we please by taking 5 sufficiently small, and the same 
being true of 2dM, it follows that | Ji | + | /a | may be made as small 
as we please by taking 5 sufficiently small and independent of x. In Ji, 
the function 

fit) -f ix) 

is evidently holomorphic for e S x ^ 1 — «, 5 ^ < ^ 1 — 5, so that con- 
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stants Ml and M2 exist such that \ <fi(x, t) \ < Mi and 

d(p(x, t) 



dt 



< M2 



for all these values of x and t. Integrating by parts, 

(p(x, t) sin (2n + l)ir(t — x)dt 

r 1 T='"* 

+ ^o^ {. 1)^ /'"' ^^%^ COS (2» + l)x(< - x)df, 



(2n + l)v. ^a 
and since | (p{x, t) cos (2n + l)7r(i — a;) 
d<p(x, t) 



we have 



af 



cos (2n + l)7r(i — x) 
2Mi + M2 



dt 

<p{x, t) I < Ml, 

< M2, 



J2 < 



(2n + l)7r • 



Consequently, after fixing a 5 so small that | /i | + | J2 1 < v, where 77 is 
as small as we please. Mi and M2, which depend on 5 but not on x, are 
fixed, and the last inequality now permits us to determine an no (dependent 
on 5 but independent of x) such that | J2 1 < »? f or n S «o and all x in 
(e, 1 — e) . The expression to the left in (37) is then less than 217 in abso- 
lute value, which proves our theorem. 

We shall now derive a formula which is useful in estimating the 
remainder in some series occurring in the next chapters: Let 

k a constant and « > but as small as we please, then 

Cn 



(38) 



Z^c^ sin {2v + k)irx 



X S 1 - €. 



sm ire 
Using the identity 

2 sin wx sin (2 v + k)irx = cos (2>' — 1 + k)Trx — cos (2j' + 1 + k)Trx, 
we have 

2 sin wx J2 c^ sin (2p + k)-KX 

n+P n+p 

= J2cy cos (21/ — 1 + ^)7r2; — 23c^ cos (2i' + 1 + A;)irx 
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n+p 

= Cn COS (2n — 1 + k)vx — X) {Cv-i — c^) cos (2>' — 1 + k)T!-x 

v=n+l 

— Cn+p COS (2n + 2p + 1 + k)wx, 

and taking the sum of the absolute values of all terms to the right and 
remembering that c^_i — c^ S 0, 

n+p n+p 

2 sin xx 2^ c^ sin (2j' + k)T!-x ^ c„ + X) (c^-i — c^) + c„+p = 2c„, 

i'=n+l v=n+i 

whence (38) follows at once. 

8. Bernoulli functions and numbers. We now apply the theorem of the 
preceding paragraph to f{x) = 5 — x; (34) gives, upon integration by- 
parts, ao = 0, a„ = 0, bn = l/»ir (« S 1), whence 

sin 2mrx 



1 



, -1^~' 0<x<l, 



the series being uniformly convergent for e ^ x ^ 1 — e. The series to 
the right converges also for x = and x = 1, the sum being zero in either 
case, while the left side becomes ± J. Let [x] be the greatest integer 
contained in x, so that ^ x — [x] < 1, then ^ + [x] — x remains un- 
changed when an integer is added to x, and equals | — xforO<x< 1. 
Since the Fourier series to the right also has the period unity, we see that 
if we write 

(39) Pi(x) = i + [x] - X, 

then for any x which is not an integer 

/ . ^s -r^ / V x^ sin 2mrx 

(40) Px(x) = Z -^^ , 

while for an integer m, Pi(m + 5) -» ^ and Pi(m — 5) -» — | as 5 -> 
through positive values. Now write 

. _ Y^ COS 2n7rx 
P2(x) - 4- 2nV2 ' 

the series being uniformly convergent for all values of x, since its general 
term is less than 1/n^ in absolute value. Then, since (40) converges uni- 
formly for e ^ X ^ 1 — e, we have 

^^^^^'^ = - P.(x) =x-h 



dx 
(theorem XXI cor.) and therefore 



P2(x) =-^-^-^+c 
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in the latter interval; to determine the constant of integration, we observe 
that, on account of the uniform convergence of the series, 

r^ -r. , s ■, r' A COS 2nwx , ^ C^ cos 2m!-x , 



and therefore 



£{^-+')'"'-o 



or c = T 



_ JL 



12- 



Now the polynomial |x^ — |x + yV takes the same value (= yj) for 
X = and x = 1; consequently the expression |(x — [x])^ — |(x — [x]) 
+ iV admits the period unity and is continuous for all values of x, so that 
the relation 

■^ ^ N , / r -.s-, 1 / r n\ . 1 '^ COS 2mrX 

P2{x) = h(x - [x]y - i(x - M) + xV = Z ^;^,;^. 

which we proved for < x < 1, subsists for all values of x since both 
members are continuous everywhere and admit the period unity. We 
observe that the integral values of x form no exception as in the case 
of Pi(x). 

Continuing this process, we define functions Pk(x), k = 2, S, ■•■, 
by the relations 

P.(x) = (-1)^'^^, 

(41) J Pk+i{x)dx = 0, 

P,+i(x + 1) = P*+i(x), 

and show by complete induction that Pk{x) is a polynomial of degree k 
in X — [x] with rational coefficients, that Pkix) is continuous for all values 

of X when /b > 1, since Pk+iil) - Pk+i(0) = (- 1)" f Pk{x)dx = 0, and 
that 

■p ( \ — Y" ^^^ 2mrx 
(42) 

P2k+l{X) = 2-1 22k^2k+1^2k+l > 

these series being absolutely and uniformly convergent for all values of x. 



ik- 
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From (42) it follows that, when m is an integer 

Pik+iim) = Pik+xim + I) = 0, 

1 °° 1 

(43) Pikim) = 22*-i^2* 2 ;^ , 

P2.(m + i) = - (^1 _ ^^P,,{m), 
the last relation resulting from the identity 

f- (-l )" of 1 f- 1 _/, 1 \f- 1 

Furthermore (42) gives 

(44) P.(l - x) = (- l)*P*(x) = Pk{- x). 

"We shall now show that /or k odd and > 1, the equation Pk{x) = has the 
roots 0, i 1 and that Pk{x) > for < x < ^, Pk{x) < for ^ < x < 1, 
while for k even, there is one and only one root interior to each of the intervals 
(0, i) and (J, 1). From (43) it is seen that Pk{0) = P*(J) = P*(l) = 
for k odd, but that for A; even none of these expressions vanishes. For 
k = 2, our theorem is true, since P2(x) = |x^ — ^x + 12 in (0, 1); we 
may therefore assume the theorem proved up to the index k — 1 and 
deduce its validity for the index k in the following manner. If, for k odd, 
Pk{x) = has a root in (0, 1) distinct from 0, |, 1, we may assume, by (44), 
that it is interior to (0, ^) . But since 

by (41), Pk-i{x) = would then, according to RoUe's theorem, have two 
roots interior to the interval (0, |), contrary to the proven facts for & — 1 
even. Thus Pk{x) does not change its sign for < x < |, and since for 
X = 0, 

^f = ^.-.<») 

is positive by (43), we find that Pk(x) > for < x < |, and (44) then 
gives Pk{x) < for ^ < X < 1. If, for k even, Pjfc(x) = has two roots 
interior to (0, |), Pfc_i(x) = would have a root interior to the same 
interval, contrary to the proven facts for A; — 1 odd. The same argument 
applies to the interval (|, 1). 

The expressions P24(m) = P^kiO) of (43) also occur in the expansion 
of the cotangent in a power series. From the decomposition of the co- 
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tangent into partial fractions we obtain 
o cot 5 = 1 - 2^ ; 



2 ""'2 „ti47rW-<2' 

the series being uniformly convergent for | < | ^ 2ir — e. For such values 
of t, we have 

and by Weierstrass's theorem on series of analytic functions, we obtain 

f f oo « 1 

i pot - - 1 —Yt^''-Y -' 

We now define the A;th Bernoulli number Bk by* 

and obtain ^ 

^Aa^ ^ J 1 ^B^t^" 

(4^) 2^«*2 = 1-S"(2ifcy!- 



In 



t_ .t__U e''/2 + e-''i^ _U fi__ 

2 ^*^* 2 ~ 2 e«/2 _ g-«/2 - 2 + e" - 1 



we now replace thy — ti; (46) then becomes 

(47) 1 - ^, _ 1 - 2 +^4^ (2;^)! > | < | < 2. 

Multiplying both sides by 

M /J; 

e' - 1 = y — 

we find 



^A;! ^- V . A;!A2'*'ri (2A;)! ) 



and comparing coefficients of i^*+^ on both sides 



1 Br ^ B2 



(2A; + 1)! 2(2A;)! (2A; - 1)!2! ' (2A; - 3)!4! 



Bz_, .i-D'B, 

rcMci ~r • • • ~r 



i2k - 5)!6! ' ' 1!(2A;)! ' 

* These numbers appear for the first time in Jacob Bernoulli's Ars conjectandi, Basel, 1713, 
p. 97, in connection with the problem of expressing the sum 1* + 2* + 3* + • • • + (n — 1)* 
as a polynomial in n. 

' Euler, L., Introductio in analysis infinitorum, Lausanne, 1748, p. 161. 
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this recurrent formula allows us to calculate successively Bi, B^, • ■ • 
which are clearly all rational numbers. The first ten Bernoulli numbers 
are thus found to be 

D_l D_l D_l R_l R 5 R — 691 

^1 — 6) -D2 — sTSj ^3 — iJ, £>4 — SO, -Dg — 66) -"6 — 2Y87) 

R— 7 R_36 1_7 R —43861 R — lJ4_etl 

■D? — 6) -D8 — -610' ) -Da T9?^ ) -''10 — ¥3 7 • 

The comparison of (43) and (45) gives, since Pk{rn) = Ph(fi) for k > \, 

P2*+l(0) = P^k+lih) = 0, 

(48) 1 ^ 1 / 1 \ 

The left side of (47) is the special case x = of 1 - [(fe'^)/(e' - 1)], 
and the latter expression may evidently be expanded in a power series in t, 
the coefficients being polynomials in x, the series converging uniformly for 
\t\ g 2x — € and | a; | < a, where e is as small and a as large as we 
please. We write this expansion 

(49) 1 - J^ = Z (- l)*('^-«/V*(x)«*. 
Since 

we obtain from (47) and (49) 

ti- i)-»/v.(.)^^ = 1 + (l:f )(- 1 +i+±i- Dfjf,), 

and upon multiplication and comparison of coefficients of like powers of t, 

(- iy<p2k{x) - - j-^ky^ -r 2(2)5; - 1)! ~ (2F- 2)12! "^ (2A; - 4)!4! 

Bzx'"-' , , , ,.k-^ ^__B'^=^=t^ 

~ (2fc - 6) !6! "*'■■■"'' ^ ^ 2!(2fc-2)! 

+ ^~ -^^'(2^)"!' 
^2k+i ^2k ^ Bix^-i ^ ggx^*^ -^ 

(- lY<P2k+i{x) = - (2A; -:|_ lyi + 2(2fcT! ~ W- lT!2! (2^- 3)!4! 

~(2r-5)!6l"^ ^^ ^ 3!(2A;-2)! 

+ ^ ^' {2k) ! 
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for A; > 0, while <pi{x) = 5 — x, so that for < a; < 1, (pi{x) = Pi(x). 
We shall now prove that for A; > 1 and ^ x ^ 1 

(Pk{x) = Pk{x). 

Since for ^ x ^ 1, the Pi(x) are completely determined by the first two 
relations (41), it is evidently sufficient to show that the (pk{x) satisfy the 
same relations. Differentiating (49) in respect to x, we obtain, since the 
differentiation of the right hand series term by term is legitimate according 
to Weierstrass's theorem on series of analytic functions, 

and comparing coefficients of i*+i 

Integrating (49) from x = to x = 1, we find 

|:(- ly^^-'-^iH^ j\,{x)dx = [^ -g7^]]~]= 0. 
whence 

<Pk(x)dx = 0. 






Outside of (0, 1), <pk(x) and Pk{x) do not coincide, since Pk{x) has the 
period unity, but <pk (x)has not. To find the difference (pk(x + 1) — (pk(x), 
change x into x + 1 in (49) and subtract: 

f;(_ !)*(*-« /2(^,(a; + 1) _ ^.(x))^ = - te'- = -Z, ""*'''* 



whence 

<Pk(x + 1) - <pk{x) = - (- !)*<*-" ''(^^)-r 

Changing k into k + 1 and adding these equations for x = 0, 1, 2, ■ ■ •, 
« — 1, we find 

p + 2* + 3* + • • • + (n - 1)* = (- l)*(*-»'2A;!(^,+i(n) - ^*+i(0)), 

which is essentially * the formula established by Bernoulli. 

' The Bernoulli polynomials (as distinct from the periodic Bernoulli functions Ph(x)) are 
commonly defined by 

0tx J^ 00 Ik 

and by comparison to (49), we readily find the relation vk(x) = (— l)*(*""/2(pt(0) — <pt(x)) 
between this definition and that of the tfext; the latter has been chosen in view of its application 
to Euler's summation formula in the next paragraph. For further information on Bernoulli 
numbers and polynomials, see Saalschiitz, L., Vorlesungen uber die BernouUischen Zahlen, 
Berlin, 1893. 
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9. Euler's summation formula. Let F{s) be a real or complex function 
of the real variable s, having continuous derivatives up to the order 2m 
inclusive, or we may suppose s complex and F{s) an analytic function; 
furthermore let S a < 1 and b > 0. In the identity 

F{s + nb) - F{s + ab + vb) =b f F'{s + tb)dt 

*Ja-i-v 

make v = 0, 1, 2, • • -, n — 1 and add, obtaining 

nF{s + nb) -JlF{s + ab + vb) = 6 E f" F'{s + tb)dt. 
Decomposing each integral by the formula 

r F'{s + tb)dt = r^' F'{s + tb)dt + Z r*' F'{s + tb)dt 
we obtain 

m-l 

nF{s + nb) -T,F{s + ab + vb) 



v=0 



n— 1 /*v+l n—1 n— 1 /»M+1 

= & E F'{s + tb)dt + & E E F'{s + tb)dt; 

in the double sum, each integral between limits p and m + 1 occurs ju times, 
namely for v = 0, 1, 2, • • •, /i — 1, so that 

Z Z I i^'(s + <&)d< = E M i^'(s + i6)di, 

or adding to the last sum the zero term corresponding to ^ = 0, and writing 
V instead of n 

n-l 

nF{s + nb) -Y.F{s + ab-{- vb) 

>.=o 

= & E ( r* F'{s + i&)di + J/ r^ F'(s + <&)d<) . 

In the integrals with lower limit v + a, we may reduce this Umit to v by 
observing that for j* S < < j* + o, [i — o] = »' — 1, or [< — a] — v + 1 = 0, 
while for v + a S < < »/ + 1, [i — a] = v, or [i — a] — i* + 1 = 1; 
consequently 

V ([i - a] - J- + \)F\s + <6)di = f F'{s + i6)d<, 
and substituting in the previous equation 
nF{s + nb) - Y.F{b + a6 + vb) = bE P ([< - a] + l)iF"(s + <6)df 



= 6 {\\t - a] + l)F'(s + i&)di. 

Jo 
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But from (39), [t — a] + 1 = t — a + ^ + Pi{t — a), and integrating by- 
parts, we find 

b Cit- a + i)F'{s + th)dt 

Jo 

= (« - o + ^)F{s + nb) - (i - a)F{s) - f F{s + tb)dt; 

Jo 

introducing all this in the preceding formula, we obtain 

Si^(s + a6 + vb) = fpis + tb)dt + (^ - a){F{s) - F{s + nb)) 
(50) "=" -^^ 

-b j Pi{t- a)F'{s + tb)dt. 

Integrating the last integral by parts, using the first and third of (41) 
and the second part of (44), we see at once by complete induction that 

"'LF{s + ab + pb) 

= fpis + tb)dt + (J - a){F{s) - F{s + nb)) 

Jo 

^ ^ + E (- l)»^*+i>«i+i6*Pi+i(a)(F('=>(s + nb) - F»>(s)) 

+ (- l)'»&2m r p^^(^ _ a)ii'(2'")(s + tb)dt, 

Jo 
S a < 1, b > 0. 

In the special case a = 0, the values of Pk+i(a) are given by (48), and 
adding F(s + nb) to both sides in (51), we obtain Euler's summation 
formula' 

i^Fis + vb) = r F{s + tb)dt + h{F{s + nb) + F{s)) 

i'=0 Jo 

(52) + E '-",0^", — fe2i-i(2?'(2*-i) (s j^ nb) - i^(2*-») (s)) 

+ (_ l)'»62'» rP2„(0/^<'"'Ks + <?')d<. 

' Euler, L., Methodus generalis summandi progressiones. Comment. Acad. Petrop., vol. 6 
(1732-33, published 1738), pp. 68-97, and independently, Maclaurin, C, A treatise on fluxions. 
Edinburgh, 1742, art. 352 and 828-30, 847. Neither author gives the remainder term. The 
proof in the text was given, in the special case a = 0, by Wirtinger, W., Einige Anwendungen der 
Euler-Maclaurinschen Summenformel, insbesondere auf eine Aufgabe von Abel. Acta Math., 
vol. 26 (1902), pp. 255-271. 
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In view of the application of (51) to the approximate evaluation of 
the sum to the left, it will be necessary to investigate the remainder term, 
or the last term to the right. The term corresponding to /c = 2m — 1 
in the sum to the right in (51) is 

(- l)'»-i&2'"-ip2„(a)(i?'(2".-i)(s + nb) - i^(2m-i)(s)) 

= (- l)^-W"P,„{a) rV(2'»)(s + tb)dt, 

Jo 

so that (51) may be written in the form 

T,F{s + ab + vb) = J'pis + tb)dt + (| - a){F{s) - F{s + nb)) 

(53) 2m-2 

+ Z (- l)»(*+»/^^+i&*P*+i(a)(i^w(s + nb) - FW(s)) + 222m, 

(54) 222. = (- 1)'"62'" r {P2,.{t -a) - P2„{a))F^'-'->{s + tb)dt. 

Jo 

Now suppose s and F{s) real, and assume that the sign of F^''"''>{s + tb) does 
not change for S i S n, so that, writing a = 1 or — 1, aF^"'{s + tb) SO 
in the interval of integration in (54). Since Pimit) takes its maximum 
value for i = 0, we have 

I P2^{t - a) - P,„{a) I ^ I PUt -a)\ + \ P,^(a) \ < 2P,^(0), 

and by (14) and (13), 

I r{P2n,{t -a) - PUa))-aF^'"^\s + tb)dt 

S r 2P2™(0)-aiF'«"')(s + tb)dt 

= 2P2m(^^^^(2m-l)(s ^ ^b) - F'-^^'^^s)), 

so that * 

(55) I Es™ 1 < 2&2'»-ip2„(0) I F(2"'-i)(s + nb) - F^'"'-^^{s) \ . 

In the two most important special cases a = and a = J, we shall now 
obtain more accurate approximate values of R^m by introducing stronger 
assumptions regarding /^'^""^(s).* Let a = ± 1, and assume first, that 

' For a = 0, a slightly more accurate formula was given by Poisson, S. D., " Sur le calcul 
num^rique des Int^grales definies," M^m. Ac. Sc. Paris, vol. 6 (1827), pp. 571-602. Other inves- 
tigations of the remainder term (in the special case o = 0) by Jacobi, C. G. J., "De uso legitime 
formulae summatoriae Maolauriniana)," Journ. f . Math., vol. 12 (1834), pp. 263-272, and Malmsten, 
C. J., "Sur la formule, etc.," Acta Math., vol. 5 (1884), pp. 1-46 (corrected reprint from Journ. f_ 
Math., vol. 35 (1847), pp. 55-82). 

» For = 0, this method is due to Sonin, N., "Sur les termes complementaires de la formule 
sommatoire d'Euler et de celle de Stirling," Ann. de I'Ec. Norm., ser. 3, vol. 6 (1889), pp. 257-262. 
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Q,2r(2m)(s -|- ii)'^ is positive and decreasing as t increases from to n + |, and 
second, that 

« E {F^^^^is + vb + tb) + jp'(2'»>(s + vb + b- to)) 

is positive and decreasing as t increases from to 5. Then, for a = 0, 

(-l)m-iii^^ = 62'"-iP2„(0)(F(2"'-"(s + nb + hb) - i;'(2'"-i)(s + hb)), 
(56) 

<h<h. 
Since a = 0, we may write (54) in the form 

(_ l)—i6-2».«22^„ = r(P2„(0) - P,M))aF^'"'^{s + tb)dt, 

Jo 

and since P^mit) has the period 1, 

(- l)"'-i&-^'»aP2« = C (P2n.(0) - P2«(0)«2i^^''"'(s + V& + to)d<. 

Decomposing the integral in two with limits 0, \ and |, 1 respectively, and 
introducing 1 — i as integration variable in the latter, we find by the 
aid of (44) that 

( _ l)".-16-2".aP2m = r (P2™(0) - P2„(0)« 2 {F (^-"^ (S + Z'b + to) 

(57) '^o v=o 

+ p(2™)(s + vb + b - tb)dt. 

To obtain an upper bound for this expression, we shall use a theorem due 
to Tchebychef: When u{t) and v{t) are continuous for a ^ t ^ b, and both 
these functions increase {or both decrease) as t increases, then 



but 



{b - a) j u{t)v{t)dt > I u{t)dt ■ I v{t)dt, 
i}) - a) \ u{t)v{t)dt < I u{t)dt • I v(,t)dt 



when one of the functions increases and the other decreases.^" Writing 

'» The following simple proof is due to Franklin, F., "Proof of a Theorem of Tchebycheff's on 
Definite Integrals," Am. Journ., vol. 7, pp. 377-379: Consider the double integral 

r r {u(t) - u{x)){vit) - v{x))dtdx = f dt f [u(t)v(t) - uit)v(x) - u{x)v(t) + u(x)v(x)]dx 

= (6 — o) I u{i)v{i)dt — I u(i)di- I v{x)dx — I v(t)dt- I u{x)dx + (6 — o) / u(x)v(x)dx, 

or replacing x by i in the single integrals to the right 

IX X ^"^'^ ~ wW)(«'W - v{x))dtdx = (6 - a) J^ u{t)v{i)di - f^ u{t)di ■ f^ v(t)dt. 

When u(t) and v(t) vary in the same sense, u{t) — u(x) and v(i) — v{x) have the same sign, so 
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u{t) = P2„(0) - Pi^it), we have, by (41), u'{t) = P2„_i(0 > for 
< i < §, so that u{t) is positive and increases with t in (0, 5), and by 
(41) and (43), 

£u(t)dt = [P2™(0) < - P2^+l{t)]'=l = ^P2„(0). 

The second factor in (57), 

v{t) = a Z {F^'"'\s + ph + tb) + F'-''"'>{s + pb + b- tb)), 

■-=0 

is positive and decreasing by hypothesis, and 

consequently, by Tchebychef's theorem, 

(58) P2„(0) raP(2m)(s + ^j)£^f > (- l)"'-i&-2"'«i?j„. 

Jo 

To obtain a lower bound for (57), we observe that u{t) is positive and 
v{t) decreasing in the interval of integration, whence 

Cu{t)v{t)dt > v{\) Cu{t)dt = P2™(0) i:«P^'"''(s + vb + P). 

But aPf^m) (g ^ j,5 ^ 1 ji, ^ ^5) ig positive and decreasing when t increases 
from to 1, and consequently 

aPem)(s ^ ^j», ^ ij»,) > r ^^^(^"^(s + ^6 + p + tb)dt, 

Jo 

so that, by the preceding inequality, 

(59) (-l)'»-i6-2'»ai?2„ > P2„(0) raP(2m)(s + l?, + ib)di. 

Since by the first hypothesis on P(^'"^(s), the integral 

f aP(2m)(s + ^5 + f5)df 

is a continuous and decreasing function oi h as h increases from to J, 
the comparison of (58) and (59) shows that there exists an h interior to 

that their product is positive, and therefore also the double integral, which proves the first part of 
the theorem. The second part follows from the first upon replacing u{t) by — u{t). 
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(0, 5) such that 

(- l)"'-i&-2"'a222„ = P2„(0) TaF^-Ks + hb + tb)dt, 

and evaluating the integral, we obtain (56). 

Now let a = ± 1, and assume first, that aF^^"''>{s + tb) is positive and 
decreasing as t increases from to » + i, and second, that 

a Z (F(2m)(s 4. p + ^?, + tb) + F(2'»)(s + hb + i'b - tb)) 
is positive and increasing as t increases from to J. Then, for a = J, 

<"»' 0<»<J. 

We obtain from (54), in the case a = i, 

(- 1)"'&-2'»E2« = r \P2m(<) - P2n.(|))-ai^(^'"Ks + ^6 + <&)df 

= r {PUt) - P2mmaZF^^'-\s ^hh + vb + tb)dt 

and decomposing the integral in two with limits — |, and 0, | respec- 
tively, and introducing — t as integration variable in the first, (44) gives 

(- l)'»6-^'»ie2. = (\pUt) - PUh))-cc"'ZiF^'"'Ks + P + V6 + tb) 

+ F(2"'>(s + 16 + j'f) - <fe))d<. 

Here the first factor under the integral sign is seen at once to be positive 
and decreasing, while the second factor is positive and increasing by 
hypothesis, and (60) is now readily proved in the same way as before. 

CHAPTER II. 

Definite integrals for the Gamma function derived from its definition as an 
infinite product. 

10. Expression of the Gamma function by Euler's integral of the second kind. 
Integrals for P{s) and Q{s). Assume 9t(s) > (whens = x + yi, 9t(s) = x 
is the real part of s), let t be a real variable, and define i' as e''°«', where 
log t denotes the real value of the logarithm of t; since 



dt\s J' 
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it follows that 



f t'-'dt = lim [-1' ' = -• 



Replacing 8 by s + 1 and subtracting the new integral from the previous 
one, we find 

and, generally, from 

'1 _ _ („ _ 1) ! 

/„ *' '^^ ~ *^" '^^ = s(s + l)---(s + n-T) 
it follows that 

j/ ^1 i)<^i- s{s + l)---{s + n-l) (s + l)(s + 2)---(s + «) 

n! 






s(s + l)---(s + «)' 
or finally, replacing i by tjn, 

n\ n' 



r'"('-s)"*= 



s(s + 1) • • • (s + ») 



For » -^ 00, the expression to the right approaches the limit r(s), as is 
seen from (J4") upon replacing n by « + 1 and a by unity, and since 
(1 — tjnY -> e~', the above equation suggests that 

(61) r(s) = r e-H'-Ht, 3i(s) > 0." 

Jo 

To prove this, we must first establish the convergence of the infinite 
integral, which is known as Euler's integral of the second kind. Writing 
s = X -{- yi, and assuming A S x S e > 0, where A is as large and e as 
small as we please, we have \e~H'~^\ = e~'V^^, and for < < s 1, 
e~ H''~^ < f '"S while for t sufficiently large, e * > f^+^ and 

Hence by theorem XVI corollary and theorem XIX, the integral in (61) 
exists uniformly for A s 9J(s) s e > 0. In the second place, it must be 
shown that the left and right members of (61) are equal, or that 

lim f t'-^(l --Xdt = r e-H'-^dt, $R(s) > 0. 



" For 5 real and positive, this expression is due to Euler (J^^). 
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We have 

e-'f-^dt - f»-M 1 -- ) dt 
Jo Jo \ nj 

and the second integral to the right approaches zero as « -> =o , on account 
of the convergence of the integral in (61). Using the inequality'^ 



<e-' - f 1 --V^ 2-e-' for ^ t ^ n, 
\ nj n 

we find 

ir('"'-('-s)")'-""hrO"'-(i-s)")""'* 



g —e-'t^'dt<- I e-H^+^dt-^O, 
Jo n nJo 



the last integral being convergent and independent of n. Thus the proof 

of (61) is complete. 

Replacing t by at, where a > 0, we obtain from (61) 

(62) ^ = r e-^H^^dt, a>0, 9?(s) > 0. 

The integral 

I e-'f-Hog tdt 

Jo 

converges uniformly for A S a; S e > 0, since for t sufficiently small, 
I e-'f-i log f I < r-i I log < I < i("2)-i, and for t sufficiently large, 

I e-'<»-i log < I < <-2 log t < r"2; 
" Proof: for 6 g a g 0, we have 

g 6» - a" = (6 - o)(6»-' + 6»-2o + • • • + a""') ^ n6"-'(b - a). 
When 6 = e~'/", where ^ / ^ », and a = 1 — (</n), we have 

so that, the series being alternating with numerically decreasing terms, 
Consequently 



— •e-'-^-e-'-^2e-'-. 
2 « 2 n re 



The inequality and its use in proving (61) are due to Watson, G. N., "An inequality associated 
with the Gamma function," Messenger of Math., vol. 45 (1915), pp. 28-30. 



78 T. H. GRONWALL. 

consequently by theorem XXI, we may differentiate (61) in respect to 
X under the integral sign, and since 

ar(s) ^ dr(s) 

dx ds ' 

we find 

(63) ^^ = r e-H^^ log tdt, 9?(s) > 0, 
and the repeated differentiation under the integral sign gives 

(64) ^^ = £ e-'i-'(log tydt, 3J(s) > 0. 

Integral expressions for the functions P{s) and Q{s) may be obtained in 
the following manner: the expansion 



(-1 )- 






converges uniformly for < e ^ i S 1, and therefore, by theorem XXI 
corollary 



f 



s + V 
For dt(s) = X > 0, the integral I e~H'~^dt converges, and since 

Jo 
[s + j-l = \x + v + yi\ ^ X + V ^ X, 



we have 



^ (- 1)" e 



V S+V 



= — ^—,-^0 as «-»0, 



^ A s + V 
so that, by the definition of P{s) (/33), 

(65) ^(«) = S ^V) = f '-''^'''' ^(^) > «• 

Since Q(s) = r(s) - P(s), (61) and (65) give 
(66) Q{s) = J" e-H'-'dt. 

As the integrand does not become infinite at the lower limit for any 
finite value of s, the last integral converges uniformly for — A ^ 9?(s) S A, 
where A is as large as we please, and consequently the entire function 
Q{s) equals the integral (66) for all finite values of s, by theorem XXIII. 
The integrals (65) and (66) are due to Legendre and Prym (J'^' ^^). 

11. The hypergeometric series F{a, /3, y; s) and its value for s = 1. 
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Etiler's integral of the first kind. The name hypergeometric was applied by 
Wallis 1' to the power series in s 

a-B a(a +l)-/3(/3+l) 
^"^1.7*^ l-2.7(7 + l) 

a(a +!)(« + 2) •|3(/3+ l)(/3 + 2) 
"^ 1.2-3-7(7 + l)(7 + 2) '^ '"' 

where a, and 7 may have any real or complex values, except that 7 
must not equal zero or a negative integer (for when 7 = — n, the n + 1st 
and following coefficients become infinite) . Using the notations 

(«)„ = «(« + l)(a + 2)---(a + n -1), 
(67) 

(a, «) = (to — l)!TO"/(a)n, (a, «) -» r(a) as n -> 00, 

already introduced in J § 14 and J § 2 respectively, we may write this; 
series 

(68) ^(«,^,,.,)=|:W|),,, 

Many of the elementary functions are special cases of this series; for 
instance 

(1 + s)" = Fi-n, /3, ^; s), log (1 + s) = sFil, 1, 2; - s), 
e' = lim^l, |3, 1; s/fi); 

the essential properties of the hypergeometric series, and its intimate con- 
nection with the Gamma function, were discovered by Gauss (J°). By 
(67), the absolute value of the general term in (68) may be written 

1(1, ^) (7,^) s- (1, '')(t, »') 



and since the first factor approaches 

r(7) 



„l+!)i(Y-a-^)> 



r(«)r(/3) 

as J- -> 00, and S l/j'i+* converges for A; > 0, but diverges for k ^ 0, it 
is evident that the series (68) converge absolutely for | s | < 1, but is 
divergent for \ s\ > 1, and that for j s | = 1, the series is absolutely 
convergent when 5J(7 — a — /3) > 0." To find the value of F(a, /3, 7; 1), 
which is a convergent series when 9? (7 — a — /3) > 0, we begin by 

'' Wallis, J., Arithmetica infinitorum, Oxford, 1656. 

" By more elaborate methods, it may be shown that for |s| =1, the series is (conditionally) 
convergent for g B{-f — a — fi) > — 1, divergent for B{y — a — /3) ^ — 1, except for s = 1, 
where the series diverges for B(y — a — 0) ^0. 
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establishing the following relation between the three hypergeometric 
series Fia, /3, y; s), F(a, ;3, 7 - 1; s) and F{a, /3, 7 + 1; s) for \s\ < 1: 
(69) 7[7 - 1 - (27 - « - ^ - l)s]F{a, ft 7; s) 

+ (7 - a)(7 - ^)sF(a, /3, 7 + 1; s) 

- 7(7 - 1)(1 - s)F{a, ft 7 - 1; s) = 0. 

In fact, the constant term on the left equals 7(7 — 1) — 7(7 — 1) = 0, 
and the coefficient of s", where « > 0, is seen by (68) to be 

^^^ ^l)n(7)n ^^^^ " ^ ^^l)„_x(7)n-l 

\ ( \l a\ («)«-! (/3) n-l 

+ (7 - a) (7 - ft - 



or denoting by 



(l)„_i(7 + l)n-l 

( ^\f (" )"(^)" («)n-l(ftn-l \ 

^^^ '^V(l)n(7-l)n (l)„-l(7 - l)n-J 



Un = 



(1)«(7 - l)n 

the coefficient of s" in F{a, /3, 7 — 1; s), 

7(7 - i)Un ( _. , n(l + « + /3 - 27)(7 + n - 1) 
y + n-lV "*" (a + n - l)(|3 + n - 1) 



w(7 - «)(7 - ft , , ^ IS 

+ (« + n-l)(^ + n-l)-(^ + "-l) 



n(7 + « — 2)(7 + n — 1) 



^i)\ 

(a + W- l)(;8 + «- 1) ; 

n7(7 - l)Un / _ 1 , (7 - «)(7 - ff) 



1 + 



~ 7 + n-lV ^(a + n - l)(/3 + w - 1) 

(« + /3-7 + «-l)(7 + n-l) \ _ 
+ (a + «-l)(^ + n- 1) j-"- 

In (69) we now let s tend toward unity through positive and increasing 
values; since F{a, /3, 7; 1) and F{a, |3, 7 + 1; 1) are convergent series on 
account of 9?(7 — a — |8) > 0, Abel's theorem on power series (J") 
shows that 

lim^a, ft7; s) = F{oc,P,y; 1) 

and 

lim Fia, ft 7 + 1; s) = F{a, ft 7 + 1) !)• 

00 

Furthermore (1 - s)Fia, ;8, 7 - 1 ; s) = 1 + Z (u, - w.-Os" f or | s | < 1 
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and since 



(1, »)(7 - 1, n) 



l+i:K-t^.-i)=lim«„ = lim („^„)(^^^) ^v- 



■«-3 



= 



on account oidt{y — a — 0) > 0, the theorem in question gives 
lim (1 - s)Fia, /3, 7 - 1; s) = 0. 

Thus we finally obtain from (69) 

F{a, 0, y; 1) = ^^^^^^nc^, /3, 7 + 1; 1) 

and replacing 7 by 7 + 1, 7 + 2, •••,7 + «— 1 and multiplying the 
corresponding equations 

(7 - a)„(7 - /3)„ 



F{a, /3, 7; 1) = 



(y)n{y - a - fi)n 
(7, n){y - a - 0,n) 



Fia, fi,y + n; 1) 
F{a, p,y + n; 1). 



(7 - a, n)(y - &, n) 

The first factor to the right approaches the limit r(7)r(7 — a — (3)/ 
r(7 — a)V{y — (8) as w -» 00, and for w > | 7 | , we have 

(l).(7 + n)J 



\F{a,^,y^n; 1) - 1 | = 



(«).(/3). 



^Z 



(«).(|3). 



~;^i(l)v(n - I7I )v 



< 



(1).(7 + w) 

I I /3 I ^ ( I a I + 1).(|/3| + 1) 



n - I 7 U=i (l).(n + 1 - I 7 I ). ' 

the last series converges for«> |7| + |a| + |;8| — 1 and is evidently a 
decreasing function of n, and since l/(n — I7I) ->Oasn^ 00, we have 
lim F{a, /3, 7 + w; 1) = 1. We therefore obtain the final result 

n— ^00 

(70) F{a, &, 7; 1) = r(7^-I)r(7-^) ' 9f(7 - « - /3) > 0. 
From this formula we shall now derive the relation 

(71) J'i-Kl - ty-'dt = f§^> 9t(«) > 0, 3t(/3) > 0, 

where the integral is known as Euler's integral of the first kind '^ 
Multiplying the binomial expansion, convergent for | < | < 1, 

1' J'" (for a and /3 real and positive). 
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by i°~^ and integrating from «tos(0<€<s<l), we obtain, since the 
series converges uniformly 

(a).(l - ^). 






= '- F{a, 1 - ^, a + 1 ; s) - ^ F{a, 1 - /3, a + 1 ; e). 

Now let € ^ and s -> 1; since 3J(a) > 0, 3{(/3) > 0, the integral tends 
toward the left member of (71), «" -> and on account of 

5{(« + 1 - a - (1 - ^)) = $R(^) > 0, 

F{a, 1 — /3, a + 1; 1) is convergent, so that Abel's theorem on power 
series gives F(a, 1 — jS, a + 1; s) -^ F(a, 1 — /S, a + 1; 1). Applying 
(70), we therefore obtain 

£t'-Ki - ty-^dt = iF(a, 1 - ^, « + 1; 1) 4^{g(^ 

_ r(a)r(ff) 
r(a + ^) • 

We may now express the hypergeometric series as a definite integral. 
Let s have a fixed value such that | s | < 1, then \ st\ g | s | for ^ < S 1, 
and the binomial expansion 

converges uniformly in t for OSf^l. Assuming 9J(a)>0, dt{y — a)>0, 
multiplying by i°~i(l — ty-"'^ and integrating from to 1, we may 
therefore integrate term by term on the right side and obtain 



Jo v=<l \^)v Jo 



dt 



_^(S± . r(a + v)Tiy - a) 
".tSd)/ ■ T(y+v) 

by (71), or since T{a + v) = (a),r(a) by (J2), 

Ct-Kl - ty—K-i- - str^dt = ^^'''^l\l~ ""h ia, &, y; s), 

(72) Jt> ^ ^y) 

I S I < 1, $R(a) > 0, 9e(7 - a) > 0. 



THE THEORY OF THE GAMMA FUNCTION. 83 

12. Application of Euler's summation formula to log T(s + a) and \p{s + a). 
To secure the single-valuedness of the logarithms with which we shall 
deal in the following, we introduce a cut in the complex plane along the 
negative real axis (from to — oo) and subject the complex variable s 
to the restriction that it shall never cross this cut. Then every branch of 
the function log s is single-valued, and we shall define log s as that 
particular branch which is real for real and positive values of s. In par- 
ticular, we have log (s)„ = log s + log (s + 1) + • • • + log (s + n — 1). 
We shall now apply (50) and (51) to the case F(s) = log s, 6 = 1 and 
obtain first a proof that (s, n) approaches a limit as n -> =o for any s 
not on the negative real axis (that is, a new proof of the convergence of 
Euler's infinite product for r(s), J3 and 3', for these values of s) and 
second, Stirling's asymptotic series for log r(s). 

Making F(s) = log s, 6 = 1, ^ a < 1 in (50), we find, since 

J"log (s + t)dt = (s -\- n)(log (s + n) - 1) - s(log s - 1) and F'is) = - , 

m-l 

log (s + a)n = JL^og {s + a + v) = {n + s + a - i) log (s + n) 

- n - (s + o - I) log s - J 'g _^ ^ dt; 

in particular, we find for s = 1, a = and replacing n by n — 1 

log (n - 1)! = Elog (1 + v) = (n - i) logn - n + 1 - J f^^df. 

Adding (s + a) log n to the last equation, and subtracting the first, we 
obtain 

log (s + a,n) = log (to — 1) ! + (s + o) log to — log {s + a)n 
/^gN = - (n + s + a - i) log (^1 -f- M + (s + a - J) log s 



+ '-X"""ff.-+X 



s + t 



Since Pi{t) = -^ by (41), an integration by parts gives 

r- Pi(f-o) ^^_P2(o) P2(a) r P2(f-o) ^^ 

Jo s + t s s + n Jo {s + ty ' 

foo p /'4 ^\ 
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verges for every s not on the negative real axis (theor. XVI cor.). Con- 
sequently 

(74) lim r?^^^^^dt= r -Pi(^-«) ^f _-P^W r^'^^~ ^^dt 

and similarly lim I ^ ', di = I ^r-r"* ^^ exists ; furthermore 

lim (n + s + a — 5) log { 1 H ) = lim {n -\- s -\- a — \). — = s, and 

substitut ng all this in (73) we find 

c" P (() r 

limlog (s+a, w) = (s+a — §) logs— s+l— I , ^ dt -{- I 



S + t 



dt. 



From the existence of the limit to the left we infer at once that 
lim {s -{■ a, n) = r(s + a) exist for every s not on the cut, and that for 

such values of s 

(75) logr(s + a) = (s + a - i) log s - s + fc' + f ^~^ dt, 

X°° P (t) 
^ '1 .dt. 

— 477 ^^ ~ I 1 \_ f dt ->0 as n -> 00 

XX p (f\ 
^ \_ . dt; making a = 0, s = n and 

s = 2n in (75), we obtain 

2 log r(n) = (2n - 1) log n - 2n + 2fc' + 2 f — x^df, 

•/o W -p C 

log r(2n) = (2n - i) log n - 2n + k' + l^ -^^^^^dt, 
and subtracting 

Jo n + t Jo 2n + t * r(2n) 

2ir(n)ni 

" ^^4(1 + 1)- ••(l + n-1) ' 

so that, as n ^ 00, fc' = log (2'^T(^)) = log V27r. This is essentially the 
method given in / §4; another is the following: making s = | ± yi, 
y > and a = 0, so that for i S 0, | s + i | > y and \ s + t \ > ^ + t, 
I s + < I 2 > 2/i(i + <)3, it follows from (74) that 



-> 00. 
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j r Pi«) P,(0) P,(0) n dt . 

I Jo i ± 2/* + i 2/ y Jo ih + ty 

Now, by (J6) 

r(i + yi)T{^ - yi) = ^ ^''' 



sin 7r(i + 2/i) e'^tt+i"' - e-'tti+y" 

^ 27r 2x6-'" 

~ e-" + e'"' ~ 1 + e-^"" 

and adding the two equations (75) iov a = 0, s = \ -\- yi and s = | — yi, 

log 2ir - Try - log (1 + e"'") 

= 2/i log (I + yi) - (^ + 2/i) - yi log (^ - t/i) - (i - i/i) 

Jo i + 2/i + < Jo i - 2/» + < 

But log (i + 2/i) = log 2/ + log i + log ( 1 + 2"^ j = log 2/ + ^ + 2^. + 
• • •, and similarly log ^ - yi) = log ?/ - 2" - 2w + " ' " ! introducing 

this in the preceding equation and letting y -> 00 , we find 2k' = log 2t. 
Thus (75) finally becomes i" 

log r(s + a) = (s + a - I) log s - s + log V2"x + a)(s, a), 
(76) pPi^-a),, 

When a = 0, 

Jo S + f „=o J^ S + < 



d< 



I^ fz-l'%rfdt = t f-i^.dt, 

y^oJo S + V + t r=oJo S + V + t ' 



and performing the integrations in the last expression, we obtain Guder- 
mann's formula for w(s, 0), (J 13). 

In (76), we now integrate the expression for co(s, a) by parts; the 
result may be read off immediately from a comparison of (50) and (51), 
ob serving th at P*(s) = (- l)*-i(A; - l)!/s*. Thus we obtain 

'« For o = and s real, this formula is due to Gilbert, Ph. (J^^), and for o = and s complex 
to Stieltjes, T. J., Sur le d^veloppement de log r(a), Journ. de Math., ser. 4, vol. 5 (1889), pp. 
425-444. For ^ a < 1 and s complex, see Lindelof, E., Le calcul des r^sidus, Paris, 1905, 
pp. 95-97. 



86 T. H. GRONWALL. 

To form an estimate of the remainder, let s = | s | e'\ — tt < 6 < tt ; 
then 

I s + < |2 = (I s I COS + 0^ + I s |2 sin2 6 

= {\s\ + ty cos^ ^+ {\s\ - ty sin^ I S (I s I + 0' cos^ | , 
and since | P2m{t — a) | S PimiO), we obtain 

(2TO - 1) \P2n,{0) r dt 



i-^^^<^-H"m^^\< "\-j',rv . 



(I s I + ty- 

^ (2m - 2) !P2»(0) 
cos^"" 0/2 • I s 1 2™-i ■ 

Joining the term with k = 2m — 1 in (77) to the remainder, we may 
therefore write 

„(, n^ - Y(- nw-)i/^i^-i)!^*+i(«) 4. (2m-2)!A.(0) h 

(78) :fe s* "^ cos^" 0/2- 1 s I 2--1 ' 

\h\ <2. 

In the particular cases where a = or a = |, P/fc+i(a) = for k + 1 odd, 
and therefore, replacing A; by 2A; — 1, (78) becomes 

(79) '^'^®' "'' ~ ^/ ^-^ s2*+i "^ cos^- 0/2 • I s 1 2™-i ' 

\h\ < 2, a = or |. 

When s is real and positive, we may use the much more accurate remainder 
terms (56) and (60). In fact, for a = and a = — 1, 

is positive and decreases as t increases from to =o , and 

a(P(2'»)(s + p + t) + P(2'»)(s + p + 1 - t)) 

(2m - 1) ! (2m - 1) ! 



(s + J- + 0"" (s + J' + 1 - ty'" 
is positive and decreases as t increases from to J (since the derivative 
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in respect to t is negative). Similarly we see that for a = 5 and a = — 1, 
the conditions for the applicability of (60) are satisfied. Taking the 
values of PikiO) and Pikih) from (48), we therefore obtain'^ 

^ -^ (- l)*-'g, 1 (- 1)"-^^^ 1 

(80) " ^"^ (2^ - 1) -2^; ■ s^"-^ "^ (2m - 1) -2771 (s + /i)2'»-i ' 

< h < ^, s > 0. 



(81) 



^ ' ""^ h\ 22*-iy (2A; - 1) -2^ s2*-i 

^V 22™-V(2m-l)-2m(s + /i)'™-^' ^^^^^> ^^"• 



The infinite series 

y;(_ ^yw^-Dm C^ ~ l)!-Pt+i(ffl) 



t=l 



the first 2m — 2 terms of which occur in (78), does not converge for any 
value of s. First suppose a + and a + 2> then the ratio of the terms 
for which A; = 2m + 2 and k = 2m may be written, using the expressions 
(42) for P2m+3{a) and P2m+i(a) : 

y, sin 2mra 
m{2m + 1) V "n^^™+' 



2x^82 -A sin2nira ' 

T »"" 

and since X) "2m — converges uniformly in respect to m for m S 1, 

1 w 

we have 23 ^ — -> sin 2ira as to -> « , so that the second factor in 

1 ^ 

our ratio approaches unity as to -> w. The first factor, however, ap- 
proaches infinity, and consequently the series cannot converge. When 
a = or a = 5, a similar argument applies to the ratio of the terms with 
k = 2to + 3 and k = 2m + 1. 

Proceeding to the logarithmic derivative of the Gamma function, 

we find, making F{s) =-,& = l,0^a<lin (50), that 

o 



1^^ = log (s+n) - log.+(i-a) (i - ^) + J"^ 



{t - a)dt 

and since the last integral converges as n -^ », and ^{s -\- a) = 

lim { log n — 23 ~x — ZT^ ) ^y {J 19'), it is seen that 

" (80) is due to Stirling (J^') and (81) to Gauss (J"), both without remainder terms. The 
form of the remainder term given in (80) is due to Sonin.* For o = 0, the remainder term in 
(79) was given by Stieltjes'* (who in this particular case obtains | A | ^ 1). 
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^(s + a) = log s — w*(s, a) 
(82) ^^ ^ _ ^^ ^1 ^ f°° Pi(f - a)df 



1 r" 

«*(s, a) = (i - a)- + 



Treating the integral in the same manner as the corresponding one for 
w(s, a), we readily obtain 



co*(s, a) = (I - a)- + £(- l)K*-'x*-2)]p(i iy.Pk(a) 

(83) 






co*(s,a) = (i-a)^+ E(-l) 
(84) 



■'■I \^ f i\\(k-nk-mii i^ ~ 1)!-P<:(<^) 



and for s real and positive 



(2m- l)!P2^(0)-/t |.| ^^ 

"^ cos^™ 0/2-1 sl^- ' l/^l ^ ^' 



(85) ^' ' 2s^ tx 2k s'"^ 2m (s + h)"-' 

<h <h 

^"^ t{\ 2^^-') 2k s'" 

(86) 

f^ 1 ^(-l)"'B„ 1 , o<h<h; 






2m (s + hy"" 

furthermore, the same argument as before shows the divergence of the 
infinite series of which the initial terms occur in (84). 
A divergent series 

«o + 7 + ^+---+^+---, 

in which the sum of the n + 1 first terms is S„(s), is said to be an asymp- 
totic expansion of a function /(s) for a given range of values oi — arc s, 
if the expansion 

Rn(s) = S-^ifis) - Sn(s)) 

satisfies the condition lim E„(s) =0 for n fixed, even though 

lim I R„{s) I = CO when s is fixed. In this sense, (78) and (79) are 

asymptotic expansions of w(s, a), and (84) of co*(s, a), for —Tr+e^O^ir — e, 
where e is as small as we please, i. e., for s inside of any infinite sector not 
containing the negative real axis. 
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It is clear from the definition that a function /(«) cannot have two 
different asymptotic expansions ao -\ + • • • and 60 + — + • • • both 

valid in the same sector (or even for a single common value of ^ = arc s), 
since then we should have 

hm s" ( ao - 60 H ^^ \- • ■ ■ -\ -^^ — I = 0, 

|sl— >» \ * * / 

whence ao = bo, ai = bi, • • • . This remark will be used in the following 
to identify various expressions for the coefficients in (78) and (84). 

13. Integrals for log r(s) and \l/{s) of Cauchy and Gauss. We begin 
by showing that 

(87) T^—dt = log- for 31(a) > 0, di(b) > 0, 

Jo t' ^ 

where log - is the principal logarithm. By direct integration we find, 

for a; S € > 0, 

f e-^^+'^^'dt = — ?-T- (1 - e-(-+'f)^), 
Jo x + iy"- 



and since 



g-(x+i!()I'| 



X + iy \ X e 

1 



-» as T -> 00 , it follows that 

(88) .['"'''^''''^'^x + iy' 

the integral being uniformly convergent for a; S e. Let dt(a) S e, 
91(6) S €, and write x + iy = a(l - u) + bu where ^u ^1, then 
X ^ e, and from the uniform convergence of our integral it follows, by 
theor. XX, that 

r --^-^-^ = r^« r (^ - a)e-(»+'^''>">'d< 

J„ a+ [b — a)u Jo Jo 

= r dt j (b - a)e-^''+'-'^'''"''>'du 

Jo Jo 



'Jo 

'» p—at 



f 

t/0 



^ ^ dt. 



But the first single integral is found by direct integration to equal one of 
the values of log-, and since the real part of the denominator in the 
integrand is not less than e for 9? (a) g e, dt(b) S «, it follows that the 
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integral is continuous for such values of a and 6. Since the integral is 

obviously real for real values of a and h, it follows that the value of logT 

in question must be the principal one, and (87) is proved. In the equation 
(J 19"), making a = 1, 

^(s) = lim (log (n + 1) - Er^) 

we now assume 3t(s) > 0; then we have, by (88) and (87), for »» = 0, 1, 
■ * ■ > 



s + 
and consequently 



~= r e-^'+-''"dt, log(n + l)= f 



»p— i -,—(«+!) I 



t 



dt, 



t/'(s) = lim ( ~ i:e-('+''>' \dt. 

Upon summation of the finite geometric series to the right, this becomes 



) = lim r ('-' - i^) (1 - e-"Od<. 



Hs. 



Now let 9'i(s) ^ € and | s | < A, where e < 1 is as small and A as large 
as we please; using the power series for e~' and e~'', we see that the first 
factor under the integral sign is expansible in a power series in t, convergent 
for I < I < r, where r is independent of s, and consequently 



-J - ^ _A <Mi = Mi(e, A) 
for ^ t < r. On the other hand, we have for t ^ r, since e < 1 



t 1 - e-' 
so that finally, f or < ^ and 3J(s) 



^e-' e-" /I 1 \ , 

- r + 1 - e-^ < I, r + 1 - e-^y ^ ' ' 



si ^^ 



< M{e, A)-e-'K 



t 1 - e-' 
Consequently the integral 

converges uniformly in respect to s for 9t(s) S e, \ s\ ^ A, and 
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as n -> 00, SO that finally 

(89) Hs) = X" (t - 1^')^*' ^^'^ > ^• 

We may transform this integral by writing f = log (1 + w , A = e* — 1, 

and since 

< J^^d< < J^y = log ^^^ -> as 5 ^ 0, 

this may be written 

or finally, replacing m by < in the last integral, 

(e-'- (1 + 0-^)7, 3J(s)>0. 
Again, replacing e~' by < in (89), we obtain 

(91) ^^'^ = £{-ikt-r^t)''' ^^^)>'- 

For s = 1, we have xl^{s) = - C, whence by (89), (90), (91) 

c-X"(r^-7)«-* 

-r(Ti-<--)f 

-.{'(r^.+i^.)^'- 

^(s) + c = l ^^^^ 

XOO y// 

((l + r^-(l + 0-')y. 3J(s)>0 

Since (89) converges uniformly for 5R(s) S €, ] s | S A, the corresponding 
integral for log T(s) may be obtained by integration according to theor. 



and consequently 

■•<»/>-' — p-" 

dt 
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XX, but it is equally simple to start from (/ 3") for a = 1, or rather the 
equivalent formula 

log r(s) = lim \{s-l) log (n + 1) - i: log ^ + '~ H . 
For 9J(s) > 0, (87) gives 
log ^; = J ^ df, log (ji + 1) = J -^ dt, 

Jr* r " H rff 

I (s - l)(e-' - e-(«+»') - E(e-" - e-('+-«0 -r 

-.'S;X[<'-'>'-+'^^]"-«->T- 

When 3J(s) S e (e < 1), I s | S A, it is seen exactly as before that 

for < ^ 0, where M depends only on e and A but not on s, and that 

the integral being uniformly convergent for 9i(s) ^ «, | s | S ^. In the 
same manner as (90) was derived from (89), we also find 

(95) log r W - f [ (, - 1).-. + "+,'4";7|'+" - ] ? . 91 W > 0.» 
14. Integrals of Raabe and Binet. The integral of Raabe" 

(96) I log r(s + a)da = s log s — s -\- log 'v'2ir, (s not on cut), 

Jo 

results most readily from (76), which gives upon in egration 

log r(s + a)da = s log s — s + log V27r + I da I ^ — j^. — . 

Jo Jo s "r t 

Here the double integral vanishes, being the limit as n -> =o of 

i«Eq. (89) and (91) are due to Gauss (J"), (90) to Dirichlet, P. G. L., Sur les int^grales 
Eul4riennes, Journ. f. Math., vol. 15 (1836), pp. 258-263 = Werke, vol. 1, pp. 271-278. The last 
of (93) was given by Legendre, A. M., Exerciees de caleul integral, vol. 2 (1817), p. 45, and (94) 
and (95) by Cauchy, A. L., Exerciees d'analyse et de physique math^matique, vol. 2 (Paris, 
1841), p. 380. 

" Raabe, J. L., Angenaherte Bestimmung der Factorenfolge, etc., Journ. f. Math., vol. 25 
(1843), pp. 146-159. 
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f r Pijt - a)dt r dt (" 

and by (41), f Pi(( - a)da = PaC* - 1) - PaCO = 0. Another proof 

Jo 

may be derived by means of (J 3') : 

r\ogT(s+a)da = lim I log r(n) + (s+a) log n— 2Zlog (s+a+i') da 
n^oo Jo L i'=0 J 

= lim log r(n) + (s + J) log n — I log (s + a)da 
»-*» L Jo J 

= lim [log r(n) + (s + i) log n — (s + n) log (s + n) 

n->oo 

+ s log s + n] 
= s log s — s 4- lim [log T{n) — (n — ^) log n + n], 

and the last limit is log V27r, as is seen from Stirling's formula. We shall 
now derive an expression for w(s, a) similar to the integral for log r(s) in 
(94). From (76) we obtain, since Pi(x) has the period unity, 

'"^^(^«^d. = iim"i:r^^(^-«) 

and from (88), for $R(s) > 0, 



co(s, a) = hm -^j-ZT-du = lim 2. „ , , . , du, 

1 r°° 

M + »- Jo 

n— 1 /*1 ••oo 

w(s, a) = lim Y. \ du \ Pi{u — a)e~'-'+''+-'>'dt. 



s + M 
so that 



By theorem XX, the order of integration in the repeated integrals may be 
inverted, and since P/x) = — ^ — xior — 1 < x < but Pi(x) = | — x 
for < X < 1, we have 



/^l />a /»1 

I Pi('M - a)e-"'du =1 (a - J - u)e-'"du + I (a + ^ - m)^-'" 

Jo Jo Ja 

as is readily seen upon integration by parts, so that 

,1 , ^ s.dt 

—at I g— (s+i-X 



du 



o(s, a) = lim Z r r(a - § - ~) (1 - e"') + e-<"] 

= lrJo h-^-.+rri=.J(i-^-"'> 



t ■ 
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Exactly as in the proof of (89), we see that for (SO, 

and that consequently, for dt{s) > 0, the limiting process may be per- 
formed under the integral sign, so that finally^" 

Writing 

(1 g-o« \ 1 



this function is holomorphic in the neighborhood of ( = 0, and conse- 

dHit a) 
quently """^r"" is bounded in this neighborhood. Using Leibniz's 

formula for obtaining this derivative, each term will contain a negative 

dHU a) 
power of ( as factor, and consequently — ^^ — -> as f -^ <» . We 

therefore have 






< Mk for t S 0, 
where Mk is independent of (. Now integrating by parts in (97) , we obtain 

co(s, a) = |_ - 1. 5(*-i TT J^^^ + ^2S^. j^ —3f^-^^ ^< 

and since 9?(s) = \ s\ cos 6 > and the partial derivatives of /((, o) are 
bounded, each term in the sum vanishes at the upper limit and 






we see that 

(99) co{s,a) = g[^.^|^^-^+ 1^,..-,^^^, , \h\<l 

Consequently (99) is an asymptotic expansion of w(s, a) for 9?(s) > 0, and 
must therefore coincide with (78), whence 

Using this in the expansion of f(t, a) by Taylor's theorem, we obtain a 
formula which is essentially identical with (49). 

*• For a = 0, this expression is due to Binet (P*). 
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It should be noted that this method of deriving the asymptotic expan- 
sion of w(s, a) only shows the expansion to be valid in the half-plane 
9?(s) > 0, whereas the method of § 12 is valid when s is not on the negative 
real axis. 

15. The generalized Schaar integral. To obtain this integral, we shall 
replace each term in the expression (/ 24) 

#(s)^|, 1 



ds ^0 (s + vY 

by an integral over a rational function (instead of the exponential used 
in § 13). Let s = x -{- yi and a; > 0; then 

r sdt ri(__i 1 \ 

X t^ + s^~ X 2\t - y ->r xi t + y - xi) 

^i r f t - y - xi t + y + xi \ 

2 Jo \it-yy + x' it + yy + x')'^'' 

= l[^^^(! + ^)^ + ^^ - * arctan*-^^- - i arctan*-* l]^^^ 

~ 2" 
Now let <T = ^ + r]i and ^ > 0, then 



whence 



2^* "'' ~ X \f + s' f + a^)'^^' 

1 2 /"" t^dt 

^ = -X (t^ + s^Xf + .^) ' 9?(«) > 0> 9f(.) > 0, 



s + 

and this is also true for o- = 0, since the integral then reduces to the pre- 
ceding one. Differentiating in respect to x, we obtain 

and this is also true for <r = 0. In fact we have ^^ 

(101) \f + s^mii'' + \s |2) cos for < g 0, s = \s\ e'\ " | < ^ < | 

and for s = I s I e'', <r = | o- | e*', -^+ e ^ 6, <p ^^- e, A S|s| S8, 
\ (T \ S 5, the integrand in (100) is consequently less than 



" Since \t'+s^\^ = \f' + \s \V«' |2 = (j2 + | 5 [2 cos 2ey + (| s p sin 2$y = {I' + \s \^y cos^ $ 
+ {t^ - \s p)2 sin^ e. 
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2At^ 
sin^ e{f + Sy ' 

so that the integral (100) is uniformly convergent in respect to s in the 
region considered, and theorem XXI is applicable; a similar proof applies 
to the case <r = 0. In order to obtain a formula applying to all values of 
s except those on the negative real axis, let us assume —ir+2e^6^'n- — 2e; 

IT X 

then we may determine an angle a such that — ^H-e^a^^ — c and 

TT TT 

— K+« = — a = o — «• In (100) now replace s by se~" and a by 

ve~", where »< is a positive integer or zero, and form the sum from v = 
to J' = n — 1 ; we obtain 






,dt. 



'=, (s + >')2 ~ TT Jo (fV" + s2)%tS<V" + v" 

But by (101), I (V" -V v^\ ^ (f^ + v") cos a'^v' sin e, and for | s | ^ 3, 

I tH^' + s^ I = I <» + s^e-'"' I s (<2 + I s V) cos {B - a) 
(102) 

S (f2 + 52) sm e; 

consequently 

I r" 2se^'<^ ^ 1 I ^ 2i s I ^ J. r" <^d< 

I j„ (f^e^' + s^)^ ,t: <^e« ^v''^^\- sin^ 6 .t:; ^-^ ' j, (i^ + 5^)^ ~* " 

as « -> 00, since the last integral converges and S — is a convergent series; 
therefore 

\at. 



V (S + J*)^ TT Jo 



(s + »')2 TT Jo (i^e'«' + s2)2;t'o<V'"' + »<' 



But from the partial fraction expansion of the cotangent it follows (see § 8) 
that 

V 2< _ 27r , 1 , 



and changing t into ie"*, we find 

2 n se'^'dt ^ r°° se'-^ndt f" se'-'H dt 

"f"^^' - TT X it'e'-' + s'Y + "^ Jo (th'-' + s'y + * Jo ith'-' + s'y e""^' - 1 ' 

The first of these integrals is found, by (100) for o- = 0, to equal l/2s^ 
the value 1/s of the second is obtained by direct integration, and the third 
converges uniformly for A S | s | g S, as is seen from (100) and the fact 
that 
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(103) 



t 



,2,r<e«» 



-1 



< 



sm € 
1 



t' sm' € 



for < i < 1, 



for t > 1.2 



Hence the application of theorem XXI to the third integral gives 



dx X fe^'^' + s^ e'"'""* - 1 ' 



tdt 



s ' 2s2 
and integration in respect to x, 

Hs) = logs - 2-^ + c{y) - 2 j" ^^ ^,-^,„,,. _ J, 

where we have denoted the integration constant by c(y). Now in 
s = X + yilet X -^ oo; then 6 = arctan y/x -> 0, so that for any a where 
I a I S 7r/2 — € and a has the same sign as y, we have \6 — a \ S ■jr/2 — e 
for X sufficiently large. Then, by (102) and (103) 



J'»00 
. 



fe^* + §2 e^""^' 



tdt 

e"-" — 1 



< 



1 (f^'^L+r ^3-^)^0 

p sm € \ Jo sin € Ji r sin' « / 



as a; -> 00 , and since also \p{s) — log s -* by (84), it follows that c{y) = 
for all values of y. Consequently 



(104) 



i{s) = logs-^-2j^ 



tdt 



<-e'" + s^ gs-'"" _ 1 ' 



< T, 



< 



- a\ < 



and upon integration by parts 



(105) '/'(s) = logs "2^ + ^ J 



1 . 1 ne'^^is'- - iV"0 



(i2g2ai ^ s2)2 



log (1 - e-^^'^Vi. 



Since \a\ < 7r/2, we have | e-^"'""* | = e-^'"^""" < 1 for i > 0, and con- 
sequently 



(106) 



log (1 - e-2-'^°*) 1 = 



» g-2ir(e" 



00 /> — 2wicOSa 



sE' 



= - log (1 - e-z-fcosa). 



the last two expressions show that log (1 — e'^"'"") becomes infinite as 
log i for t -^0 and vanishes as g-S'icosa ^j^gjj ^ _^ oo. Together with 

22 Since | e'^'"^' - 1 I g 1 e'^'""' | - 1 = 6="' «« « _ i > g2»( sin . _ j > either of 2Tri sin e or 
(27rt sin e)^ 
3! 
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(102), this shows that the integral in (105) converges uniformly in respect 
to s in the region previously considered, and we may therefore integrate 
(105) in respect to x, obtaining 

log r(s) = (s - §) log s - s + c{y) - -| ^^M-^. log (1 - e-^'"-')dt 

From (102) and the preceding remarks on the logarithm it follows that 
the last integral approaches ^ro as x -> », and from (76) and (78) 
logr(s) — (s — 7) +s ->log 'V27r, so that c{y) equals the latter constant 
and we finally obtain 

log r(s) = (s — I) log s — s + log V2ir 

\e\<T, |al < J, |0-a| <|, 

where, by (76), the last integral equals w(s, 0).^' 

Making s = 1 and a = in (104) and (105), we obtain two expressions^^ 
for Euler's constant 

^ ' ^ Jo ^^H-le^^'-l 
(108) _ ^ 

To obtain the asymptotic expansion of a)(s, 0) from (107), we expand the 
first factor under the integral sign in negative powers of s : 

„gat m-i ( -i\k—lMk—3)<ii^2k—2. ( J\m— lg(4m— 3)«<^2m— 2 

_ 2^ „2A-1 I 



^2g2«i ^ s2 ^^ s2A-l 1 s2m-3(^2g2.< _^ ^2) 

and obtain 

m-l ( 1 "\ ft 1 /'M . 

w(s, 0) = E -Vft^ • - e^^^-^^^V-'- log (1 - e-^-'^'Vi + B-, 

t=l S T^ Jo 

( ■\\m 1 /»oo „(4m— 3)a»/2m— 2 

From (102) and (106) it follows that 

^' For a = 0, 8{(s) > 0, (107) is due to Schaar, M6moire sur les integrales euWriennes et sur 
la convergence d'une certaine classe de series, Mem. Ac. Belgique, vol. 22 (1848), pp. 3-25. 

^< Poisson, S. D., M^moire sur les integrates definies, Journ. de I'ficole Polytechn., cahier 
18 (1813), p. 305. 
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or introducing t cos a as integration variable 

1 



Rm = 



(;Qg2m-l ^ (jQg (•^ 



T - r - <2".-2 log (1 - e-^")dt. 

— a) IT Jo 



Hence I s \^'"~^ | i2„ | -> as | s | -> oo and the expansion above is asymp- 
totic, so that it must coincide with (79) for a = 0. Expressing PikiO) 
in terms of B* and identifying coefficients, we therefore obtain 

(1^^) - {2k-i)-2k = ~ir ^'"'''^'t'"' log (1 - «-''"'V<; 

making k = m, a = and introducing in the last expression for | fim |, 
we find 

(110) I ^^ I = (2^ - 1) •2m * cos^™-! a cos (0 - a) * \s p"""! * 

By choosing a for a given 9 so as to make the second denominator a 
maximum, this form of the remainder will frequently be much smaller 
than the one in (79).^^ 

16. Sonin's form of the remainder in Stirling's formula. In the re- 
mainder term in co(s, 0) obtained at the end of the preceding paragraph, 
we now assume s real and positive and make a = 0, whence 

1 1 r* i2m-2 J 

or writing t = su, 

For s > 0, M > 0, < e-^"" < 1, and 

1 <X) p — 2nn'su 

log 1 _ e-2-.M = \ ^^ 

00 y,— 2nn-(8+A)u 
= Y\, . e2n.rA„. 

assuming X > 0, e^"'^" is positive and increases with n, whence 

1 00 g— 27Hr(8+X)« 

iogr^-^3-«>^^"^"-?— ^— 



2,X„1 1 

— e ' " log J _ g_2»(a+A)„ • 

Now assume X so large that for every w S 

" For further details on this point, compare Lindelof, Calcul des r^sidus, pp. 98-102. 
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(Ill) e^"^" ^1 +u^, 

then the last expression for R^ gives 



i r 

ttJo 



u""-' log YZI^z^udu > (- l)'"-iB„ 



>iJ\-Mogj— ^i,^„du, 



and introducing the new integration variables su in the first and {s + \)u 
in the second of these integrals 






1 /••» 1 



Making k = m and a = in (109), it follows that 



(2m - 1) •2m s^"-! ^ "■ ^ *" -^ (2m - 1) •2m (s + X)2"'-i " 

To secure as close a limitation of Rm as poss ble, we have to find the 
smallest value of X which will satisfy (111). It is clear that for this 
smallest value of X, the equality sign in (111) must obtain for one or more 
values of u besides m = 0, and consequently the smallest admissible 
value of X equals the maximum of 

This maximum must satisfy the equation 



let 
so that 



d(2^ = - ^log (1 + u^) + :r4r~2 = 0; 

,d(27rX) 2m2 ,, , „ 



(1 + u^y' 



and consequently /(«^) increases for m < 1 and decreases f or m > 1. But 
for u small, f(u^) = u^ + ■ ■ ■ > 0, so that finally f{u^) = can have 
only one root, for which u > 1. The numerical calculation of this root 
gives I + u^ = 4.92154 • • • and the corresponding value of X is 0.12808- • •, 
so that we may now replace (80) by the following: 
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7/t~l 



(112) co(s, 0) = 2. Toir-rrTToi; • ^*=i + 



rl (2A; - 1) -2^; s^"-' ^ (2m - 1) •2m (s + /i)^'"-! ' 

< A < 0.12808- ••, s > 0. 

r(s) 

17. Kummer's and Lerch's trigonometric series. The function ^,^ _ — c 
is holomorphic for < s < 1, and becomes infinite as log 1/s when s -> 

X^l r(s) 
^^S Y(-i _ \ ds exists. Conse- 
quently, this function is expansible in a Fourier series (35) uniformly con- 
vergent for € ^ s ^ 1 — €, where the coefficients are given by (34) : 

f' r(s) r^ r(s) 
^^S T7i s: • cos 2nwsds, 6„ = 2 I log :^^ r sin 2nirsds, 
I \i S) Jo 1 \i S) 

and replacing s by 1 — s in the expression for a„, it is seen that 

r^ r(i — s) 

a„ = 2 I log — :pr/"T — cos 2nwsds = — a„ 
or o„ = 0. To calculate b„, we observe that, by (J 3'), 

-i^'-^ - lim m^'- . (l-«)(2-«)---(m-s) 
r(l - s) " l'^„™ s(l + s) ■ • -(w - H- s) ' 

so that 

[/»! m /^l 

log m • I (2s — 1) sin 2nwsds + 23 I log {f — s) sin 2n7rsds 

m-l .-%! -| 

+ 12 \ log (»< -f s) sin 2mrsds . 

y=0 Jo J 

Integrating by parts, 

r 1 

I (2s — 1) sin 2nTsds = — - — , 

Jo W" 

and taking p — s and p + s respectively as new integration variables 
XI I log {y — 8) sin ^n-w^d^ = — 2 I log s sin 2mrsds 

v=l Jo f=l •/»-— 1 

= — I log s sin 2mrsds, 

Jo 

m~l />1 m— 1 /^r+l 

XI I log (v + s) sin 2n7rsds == XI 1 log s sin 2n7rsds 

= I log s sin 2mrsds, 

Jo 
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whence 

ibn = lim — ~ 2 I log s sin 2mrsds , 

m->oo L ^T Jo J 

and integrating the last integral by parts 

,, 1 ,. r f"! - cos2n7rs , , 1 

*^"=^iT-LJo ^ ds - log m\ 

1 ,. r r^'H — COS 4TOnirs , , , , 1 

= :^"^\ I i "« ~ log mn + log n , 

SO that finally, replacing mn by m in the limit, 



where 






§&„=— (logn + X), 



„ ,. r C^l"^ 1 — cos 4TOirS , , "I 

K = lim I as — log m . 

m->oo \_J(3 S J 



To calculate K, we shall replace the integral by one which may be eval- 
uated in finite terms. Observing that 

. . 1 TT 



s sin -KS 

is holomorphic for S s ^ §, we may write 
p/^ 1 - cos 4otxs _ p' 



sin xs 



Now 

" s 



£1/2 /»l/2 

<p(s)ds — I ^(s) cos AniTsds. 



J ''1/2 




log 



tan^ s 






= log § + log2 = log 2^-3 log 2, 



and integrating by parts 

Xi/2 I r -|»=i;: 

^(s) COS 4:mirsds = j — (p(s) sin 4:mirs 

Armr 



1 fi'^ 

I (p'{s) sin 4rmrsds 

Jo 



as m -> 00 , 



since ^(s) and ^'(s) are holomorphic, and therefore bounded, in the interval 
of integration. Consequently 



THE THEORY OF THE GAMMA FUNCTION. 103 

K = lim Fx r'^-^=-^:^ds - log m] + log 2x - 3 log 2, 



and since we have 

1 - cos4mirs o v- • ^o i\ 

; = 2 2^ Sin {2v — l)7rs, 

sin TTS t^i 

the integral is readily evaluated, giving 

[2m 2 "I 

Z o T - log TO + log 2ir - 3 log 2. 

By the definition of Euler's constant, 

4m o 2m 1 

Z- - 2 log 4to -> 2C, E - - log 2to -» C 

as TO -> 00 , and by subtraction, 

2m o 

2 7^ T - log TO - 3 log 2 ^ C, 

►=i •^'' — -i 
so that finally 

K = C + log 27r, 
and consequently 

n r(s) ^ logn + C + log27r . n^ /i 

I log fiT^j r = 2^ ^ — Sin 2mrs, < s < 1. 

^ *r(l — s) T n-n- ' 

But by (40) we have, for < s < 1 

E^^^^sin 2n7rs = (C + log 27r)(i - s), 

and furthermore, by (J 6), 

r(s) , , r(s)2 , „, s , 1 , sin ITS 

whence we obtain Kummer's series^* 

log r(s) + § log'-^ + (C + log 2t)(s - 1) = E^-^sin 2nxs, 

< s < 1, 
the series being uniformly convergent for « ^ s S 1 — e. 

The differentiation term by term of the series (113) gives a divergent 

" Kummer, E. E., Beitrag zur Theorie der Function r(x), Journ. f. Math., vol. 35 (1847), 
pp. 1-4. 
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series; a trigonometric expansion of i^(s) is however furnished by the 
following theorem:^' 
Let 

f(x) = X^ — sin 2'mrx 
»=] ^ 

be convergent f or < a; < 1 ; then the derivative of f{x) is given in this 
interval by 

fix) ^- ='£iicn - Cn+i) sin (2n + 1)tx, 

where Co = 0, provided that the last series converges uniformly for 
e ^ X ^ 1 — €, where e is as small as we please. 
Writing 

00 

g(x) = Z) (c„ - Cn-\-i) sin (2n + l)irx 

T!=0 

= 12(cn — Cn+i) sin (2n + l)7ra; + Rn, 

= 12cn sin (2n + l)7ra; — 23c„ sin (2n — l)wx + Rn 

= IZc„ • 2 sin Tra; cos 2mrx — c^r+i sin (2iV + 1)tx + Rn, 

where, on account of the uniform convergence, | Rn \ may be made as 
small as we please for 6 ^ a; g 1 — e by taking N sufficiently large, 
we have 

T , , ^ _ o sin {2 N + l)i^x . wRn 
~ 9W = l^Cn • 2t cos 2nirx — wCn+i ^ r -• , 

and integrating between the limits e and x 

I -. q(x)dx = 2Z — (sin 2mrx — sin 2mre) 

rxsin(2iV+J>£ , , rjrRn 
^ J, sm TTo; Je s" 



sin TTX 
Integrating by parts, 



dx. 



TT sin (2^" + VjTTX ^^ _ ^+1 f cos (2Ar + 1) t€ 

sin 7r€ 



p7r sin(2i\r+ l)7rx ^+1 r 



cos (2iV + 1)tx f'lr cj)s { 2N + l)7rx cos ttx 



-X 



sin ttx J, sin'' irx 



" Lerch, M., Sur la differentiation d'une classe de series trigonom^triques, Ann. Ec. Norm., 
ser. 3, vol. 12 (1895), pp. 351-361. 
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and here the expression in brackets is evidently bounded for e ^ x ^ 1 — e 

i all N, so that < 
Furthermore, 



and all N, so that our integral vanishes as N -^ «>, since then r)jJ'V_ -. -»0^' 



I -. — -- dx < -. I \ R^ \ dx -> as N -^ <x> 

I J, sin irx sin ire J, ' 

on account of the uniform convergence of the series g(x), and we therefore 
have, letting iV -> oo , 

^jj^-^Sr(x)dx = E^Csin 2mrx - sin 2nire) = f{x) - f{e), 

whence our proposition follows by differentiation. 
In the special case c„ = log n/w, 

1 °° n 

g{x) = - Z log z-XT • ^^^ <^2n + l)irx, 

and since log is positive and decreases as n increases, we have 

by (38) 

n + 1 



Z) log sin {2v + l)xx 



sin ire 



log 



so that the series g(x) converges uniformly for e ^ x ^ 1 — e. We may 
therefore apply Lerch's theorem to (113) and obtain Lerch's series (1. c.) 



(114) 



TT 

V'(s) sin ITS + K cos TTS + (C + log 2t) sin its 



= Z log ^ , , • sin (2n + l)7rs, < s < 1, 



the series being uniformly convergent for € ^ s ^ 1 — €. It should be 
noted that this series is not a Fourier series according to the definition 
in §7. 

CHAPTER III. 

The Gamma function defined as a definite integral. 

18. The Euler integrals of the first and second kind and their elementary 
properties. In the preceding chapter, the various integrals connected 

^ This is an immediate consequence of the convergence of the series for /(x) in (0, 1) and the 
following theorem, due to Cantor: When S (a» cos nx ■{- 6„ sin nx) converges in an interval 



(a, 6), then On— >0, 6„^0 as n— > oo. We omit the proof, since we shall apply Lerch's theorem 
only to the case c„ = log n/ir, where it is evident that Cn+i/(2n + 1) — >0 as re— > oo. 
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with the Gamma function were derived from the definition of r(s) as 
Euler's infinite product. As we shall now proceed to show, the same 
results may also be obtained by defining T(s) by means of Euler's integral 
of the second kind (61) 

(115) r(s) = e-H'-'dt, &J(s) > 0. 

Jo 

As shown in § 10, the integral converges uniformly for e ^ dt(s) ^ A, 
where e is as small and A as large as we please, but does not converge for 
9?(s) < 0. Since therefore s must be restricted to the half-plane dt(s) > 
throughout the present chapter, certain formulas (for instance, the asymp- 
totic expansions in § 12, which were there proved for all s not on the 
negative real axis) will have a larger region of validity than that shown 
by the proofs to be given in the following. Integrating by parts, we find 

J'^OO p "I ( = 00 /><» /»00 

e-H'dt = - e-'t' + s \ e-H'-^dt = s I e-'t'-^dt, 
L J<=o Jo Jo 

or 

(116) r(s + 1) = sr(s), 

which is one of the fundamental properties of the Gamma function (J4a). 
By direct integration, (61) gives r(l) = 1, whence from (116), for a 
positive integral n, T(n) = (n — 1)! The formulas (62) and (63) 

(117) ^ = r ''""^'~'^^> «>^' 

(118) ^^-^ = P e-H'-'- log tdt 

are derived from (115) exactly as in § 10. 

Euler's integral of the first kind is defined as 

r i»-i(l - ty-Ht, 9i(a) > 0, $R(;8) > 0; 

Jo 

t U 

introducing the new integration variable u = ^ _ , , whence t = ^ . , 
we find 



(119) 1^^-^(1-0^-^^^ = for 



u) 



>+s 



du. 



To express this integral in terms of the Gamma function we write 
a = l+u, s = a + P in (117) : 
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r(a + fi) 



1 e-»+« 



>'i»+*'->di = 



(1 + uy+» ' 
multiply by u''~^du and integrate from to « , which gives 

£ du £ e-»+») V-if«+^-idf = r(a + ^) £" (1 + J).+g ^«- 

If we replace the integrand to the left by its absolute value (which amounts 
to replacing a and /3 by their real parts), the new repeated integral con- 
verges, since it equals the right hand member with 9t(a) and 9i(/3) sub- 
stituted for a and /3; by theorem XXII, we may therefore reverse the 
order of integration in the repeated integral and obtain 

£dt J°°e-(i+")'u"-ii«+^-i du = r(a + /3) £ ^^ XuY^' ^^- 
But by (117) 

r , ,j r(«) 

and consequently 

T{a)J\-'t'-Ht = r(a + P)£jY^:^,du, 

or finally, using (115) 

(120) X\TT9^^^ = r(-«T^' ^(«)>^' ^(^)>'' 

which by (119), is equivalent to (71). 
To prove that 

(121) r(s)r(i - s) = -^^^ , < m{s) < i, 

V / V, / \ ^ sm ITS ' ^ ' 

(which is J 6 with a restriction on 9i(s) imposed by the method of the 
present chapter), we make a = s, /3 = 1 — s in (11 7), so that 

r(a + p) = r(i) = 1, 

and show that 

/'OO ,.8-1 _ 

(122) ir+^^- = ^r^s' 0<9i(.)<l. 

The integral to the left converges uniformly for e S 9t(s) ^ 1 — e by 
theorem XIX, since j u'-'^Hl + «) I S u'-^ for < m < 1 and | ■m»-V(1 + «) 1 
^ M~'~^ for M > 1, and the integrals 
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both converge. By theorem XXIII, our integral therefore defines a 
function of s which is holomorphic for < 3t(s) < 1, and the expression 
to the right in (121) is also holomorphic in the same region. To show 
that both sides in (121) are equal throughout this region, it is consequently 
sufficient to prove their equality for an infinity of values of s having at 
least one limiting value interior to the region, for instance the values 
s = (2p + i)l2q, where p and q are any positive integers such that p < q. 
Making u = f^ in (119), it is therefore sufficient to show that 






dt = 2i , , ,,„ dt = 



Ul+i^« Jo 1 + t'^ . 2p + l 

«''^~2^^ 

The roots oi 1 + f" = are t„ = e<2"+"'"'24 andi„ = e-<2"+i)"/29 (n = 0, 
1, ■ ■ •, q — 1), and decomposing into partial fractions, we have 



t^" _g/ Ar. An \ 

n=0 \' 'n f tn / 



1 + f^ 
where 

t ^p 1 - 1 - 

"~2g<„2«-i~ 2q" ' " 2g'" ' 

Now t„t„ = 1, 



is real, and 



hence 



2« + 1 

tn + tn = 2 cos ^--~ IT 

tn -~tn . 2ft + 1 ^ - 



An 1^ An \An "r An)t (^Antn "T -Antn) 



t-tn t-tn f" - {tn-\- Qt + 1 

1 (An + An)i2t - tn -In) 



and 



t' - (tn + tn)t + 1 

"'('-'-r-y+W")' 



1 {An - An)itn - tn) 

"^2 
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«-»«, J-t \l — In t — t„/ t-*» V + (<„ + tn)t + 1 



tn tn 



X 



2i 



t - 
arctan 



'« 'fl 



2i 



= irt(A„ - A„) = 2x9? (iA„) 
SO that finally 

^ / o28(2p+l)»,724 _ 1 \ 

= _«?(_ ,-e(2p+l)»i/28 . E t \ 

q \ g2(2p+l)«ri/29 — 1 y 

= £«>/'_• -1-1 \ 

g \ * ' e(2p+l)».72« _ g-(2p+l)ir</2« J 

IT 1 



ff . 2p + 1 ' 
sm — ^ — TT 
2g 



which proves our proposition.^' 

*' From (121) we may obtain the expansion in partial fractions of ir/sin vs. For decomposing 
the integral in two with limits 0, 1 and 1, oo, and changing u into Iju in the second of these 
integrals, we obtain 

sm xs ''o 1 4- « ' ^ ' ' 

and since 

1 n— I «#n 

= 2 (- D-w + (- l)» 



it is seen that 



sinirs 
For € ^ 3i(«) ^1 — «, we have 



n-l (■_ iy n-l (•_ ii>. />1 „<+»-! 4. «»-» 



ol+« |Jo^ J n + 1 ' 

and consequently 



^ (-l)" I 2 (-1)" ^ (-1)- 



sin ITS t ^ + " " + 1 — » -a s + V ' 
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In (120), we now make a = = s, obtaining 

r(2s) - j„ (1 + uy'^'' + X (1 + u)^-'^" - ^ j„ (1 + «)^''^^' 

and writing u = (1 — t^'^)l(l + fi'^) in the last integral, 

r(s)^_ 1 r',-..i .^^.^ 1 ra)r(s) 



= ^i£t-Ki-t)-'dt=^, 



r(2s) ~ 22'-i Jo " ^" "^ "" 22'-i r(s + 1) 

by (119) and (120). Since (121) gives r(|) = a/^, we obtain Legendre's 
formula (J") 

(123) r(s)r(s + 1) = 2^ir(2s), $r(s) > o. 

19. Integrals for ^(s). The infinite product for r(s). In (118), we 
replace log t by the expression 



log t = f 



CO y,— « 






/o u 

obtained from (87), and find 

d< I e-'<'-» du= dt I f{t,u)du; 

since e"" 5 e~'" as < S 1, we have 

I fit, u)\ = ± e-'^o-i ^"" ~ ^"" (± as < 5 1), 
so that 

Cdt flfit, u)\du = f dt fl f{t, u)\du-\- { dt ( \f{t, u) \ du 

Jo Jo Jo Jo J\ Jo 

dt I e-'<»('>-'^ —du 

Jo ^ 

+ \ dt\ e-H^^'^-^~ —du 

Jl Jo U 

and performing the integrations in respect to u 

f dt f \f{t, u)\du r e-'pc>-' log t dt + f e-H""-'^-^ log t dt. 

for < 9i(«) < 1, but since the infinite series converges uniformly in any finite region of the 
s-plane to which the points s = 0, ±1, ±2, • • • are exterior, and consequently the series is a 
function of s holomorphic in this region, the above equation is vaUd for aU values of s distinct 
from 0, ±1, ±2, • • •. Conversely, we may deduce this partial fraction expansion from that of 
the cotangent by means of the identity 2/sin vs = cot ir«/2 + cot [x(l — «)]/2, and then the 
above expansion of the integral gives iis a new proof of (121). 
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Ill 



dt I \f(,t,u) \du 

converges and by theorem XXII, we may reverse the order of integration 
in the repeated integral for T'{s), whence 



du I t'-' 

Jo 



g—U—t g— (l+uX 



U 



dt 



by (115) and (117), or finally, writing t instead of u in the integral, 

(124) ^(s) = ^=£ (e-' - (1 + 0-) y , ms) > 0, 

which is identical with (90). We now define Euler's constant by the 
relation 

(125) - C = ^(1) = r'(l) = £ (e-' - (1 + 0-0 J, 
whence by subtraction 

Hs) + c = f\{i + 1)-^ - (1 + t)-')j, 

and replacing 1 + < by e' 

(126) 4'is) + C = f' ^'iZ^I^' dt. 

Differentiating under the integral sign. 



(127) 



where evidently 



1 - e- 



< M iort^O, so that for 9i(s) S e 



r\^---f. 



M 

Mte-"dt = ^ . 

€^ 



Consequently the integral in (127) converges uniformly for dt(s) S e, 
so that the differentiation under the integral sign is legitimate by theorem 
XXI. We have 



J. — e ^=0 



te- 



1 ;^, = Eto-W)+'')', 



and since 



Jo 1 - e-' 



1 - e-« 
d< converges for 9J(s) > 0, theorem XXII cor. 
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gives 

v=0 Jo 

or evaluating the integrals by means of (117) 

where the last series is uniformly convergent for 3J(s) ^ e, | s | S A, since 

s - 1 I A + 1 A + 1 



1 + l* S + J' 



(1 + «.) I s + H (1 + ")(« + ") 

and consequently, remembering that ^f'(l) = — C, 

\A(s) + C = Zfrx i- V 

This may also be written 

diogr(. + i) g./i_ 1 y^|g./g_ s+^x 

the last series being uniformly convergent as before, as is readily seen by 
expanding the logarithm, and since for s = 0, V{s + 1) = 1, it follows that 

log r(s + 1) + Cs = |:(? -log(l +?)). 

For s = 1, we obtain, since r(2) = 1, 

which is the definition of Euler's constant given in J § 1, and passing 
from logarithms to numbers, the preceding equation gives, since 

r(s + 1) = sr(s), 

I^ = e-.«^(l+5^)e-^ TO>0. 

Except for the condition 9J(s) > 0, this is the product formula (/ 5). 

20. Integrals for log r(s). Raabe's integral. In the same way as in 
the case of (89), we may show that the integral (126) converges uniformly 
for 9i(s) ^ €, I s I = A, and consequently 
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(s - 1)^(1 + u{s - 1)) + C(s - 1) = fis- 1) 

Jo 



g—t _ g— «-u(«— 1)< 



1 - e 



,—t 



dt 



converges uniformly for 3J(s) S «, ] s ] S A, S « S 1. By theorem 
XX, the integration in respect to u between the limits and 1 may there- 
fore be performed under the integral sign to the right and gives, since 

r(i) = 1, 

log r(s) + C{s - 1) = J^ ^ I _e-t 7- 

For s = 2 we obtain, since r(2) = 1, 

and multiplying by (s — 1) and subtracting from the preceding equation, 

(128) log r(s) = X" [ (s - 1)«-' + ^7-"e-'"' ] f ' ^(') > ^' 

which is identical with (94) and is seen, as in the case of (89), to converge 
uniformly for 9t(s) S e, | s | S A. We shall now evaluate Raabe's 
integral 

J log r(s + a)da, dt(s) > 0, 

where, on account of the uniform convergence of (128), log T{s + a) is con- 
tinuous for 9t(s) >0 and OSagl. From log T(a) = log r(l + a) — log a 

and the existence of I log ada it therefore follows that I log T(a)da 

Jo Jo 

exists. The integral 

I = I log T(us + a)da, 

Jo 

where 9t(s) ^e, \s\^A, S^uSl{8 being arbitrarily small) may be 
differentiated under the integral sign, since the resulting integral 

dl 

du 



Jo 



converges uniformly for the values of s and u considered. But integrating 
in respect to a and applying (116), we find 

81 

Q- = s(log T(us -|- 1) — log T(us)) = s log us, 

whence integrating in respect to u between the limits and 1, 
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J' log r(s + a)da — I log T(a)da = | s log iisdu = s log s — s. 
Now, replacing a by 1 — a, 

I 'og T{a)da = I log r(l — a)da 

Jo Jo 

and consequently, applying (121) 
r log T{a)da = | f (log r(a) + log r(l - a))da 

Jo Jo 

= i I (log ir — log sin ira)da = log '^ir — ^ j log sin irada. 

Jo Jo 

To the last integral, we apply the formula sin ira = 2 sin -2" sin ^ > 

whence 

I log sin Toda = log 2 + I log sin -^-da + j log sin ^ da 

= log 2 + 2 I log sin rada + 21 log sin ( h + '^''^ ) ^ 

= log 2 + 2 I log sin irada + 2 I log sin irada 
Jo Ji 

= log 2 + 2 I log sin vada, 

Jo 

J log sin Trada = — log 2 

I log r(o)da = log Vx + I log 2 = log A/27r, 
Jo 

and finally 

(129) \ log r(s + a)da = s log s - s + log a^, «R(s) > 0. 

Jo 

We may obtain another expression for Raabe's integral by replacing s 

by s + a in (128) : 

9i(s) > 0, go < 1 



so that 
and 
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and integrating in respect to a from to 1 ; the integral to the right being 
uniformly convergent in respect to a from S a ^ 1, we may perform 
the integration under the integral sign and obtain 

(131) f\ogT(s + a)da=£[{s - i)e- +^' - ^-^]^*, 9t(s)>0. 

21. Binet's integral and the asymptotic expansion of log r(s+ffl) and if'(s+a). 
From (130) we subtract (131) and find by the aid of (129): 

log r(s + a) — (s log s — s + log V2ir) 
from (87) we obtain 



log s = I 

Jo 



dt, 



t 
and multiplying by a — J and subtracting from the preceding equation, 

log r(s + a) = (s + a — I) log s — s + log V27r 

Rewriting equations (97) and (98), 

(132) „(,.„) = jr-(„_i_| + ^)!:;:^, 

(133) /(,,„).(„_ I _l + j^,).|, 

we find as in § 14, upon successive integration by parts in (132), 

(134) co(s, a) = [ - L g^^x ^ J^^^ + -^^l Qf.,n-. e -dt. 

It is now proposed to show that this formula yields the asymptotic 
expansion (78) for 9i(s) > 0, and the proof will be conducted without 
any comparison of (134) to (77) (since we have not established the latter 
formula by the methods of the present chapter). To this purpose, we 
note that /(<, a) is a function of o holomorphic for ^ a S 1, and may 
therefore be expanded in a Fourier series uniformly convergent in this 
interval 

/(<, o) = -^ + 53 (a„ cos 2mra + 6„ sin 2nira) 
•^ 1 

= f + Z (A„e2»-- + A-ne-^"--), 



116 T. H. GRONWALL. 

where, by (34), 

^n = i(«n-*n) = fVc. a)e-''''"'da, A_„= |(a„ + *„) = f fit, a)e2»"^do, 

From (133) we find 

«o = 2 I f(t, a)da = 



and 



^-»=7r(«-^-i+r?e-o^^"'^"'« 



f \_2nin 



ae 



l^J-- 



]a=l ir 1 g(2niri— Oo"1<l=l 

„=o I L2n7rt - < ■ 1 - e-' J„=o 

Tt - t) ^ ~ : 



t\2niri 2mri — tJ 2nin{2nin — t) ' 

An is obtained from this by changing the sign of i, so that 

(135) fit, a) = - E [2niri{2nid + t) + 2nin(2niri -t)) ' 

Differentiating term by term in respect to t, which is legitimate since 
the series obtained are uniformly convergent for < S 0, it is seen that 

U^t>; 3^t-i - (-« i; • £j \^2n7rtX2nxi + f) * 2n«(2nxt - V 

and from the uniform convergence for < ^ it follows that — Qtk-i — ~^ 
as < -» 00 . For < = 0, we obtain 

rd*-^(^a)1 _ ,, ^ (- i)t-ie^»'i° + e-^"''° 

L a<*-i Lo" "^ ^ ^'nt; (2n«)*+i 

and comparing the cases k odd and A; even separately to (42), we find 

Since | 2md ± < | S 2nv, it follows from (136) for A; = 2m - 1 that 

^7f^-S ""^ I < (2m - 2) ! E (2^„ = (2m - 2) !P..(0), 
and for s = I s I e** (- ir/2 < 6 < ir/2 since 9e(s) > 0), 
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Qtln-2 ^ e-''dt\^< (2m - 2)!P2.(0) J^ e-'"-'*'di 

= (2m-2)! |^r^"\ . 
^ I s I cos 6 

Introducing all this in (134), we finally obtain 

u[s, a) - Z. ( Ij ^, + I s |2— 1 cos e ' 

(137) 

\h\<l, 0^a<l, -|<0<|, 

which, except for the slightly different form of the remainder term, is 
identical to (78). 
We have 

w*(s, a) = log s - i^(s + a) = (I - a) - - "^^^ °^ , 

and by theorem XXI, it is legitimate to differentiate in respect to s under 
the integral sign in (132), whence 

(138) aj*(s, a) = {^ -a)\+ C tf(t, a)e-'*dt, 9?(s) > 0. 

s Jo 

Integrating by parts, and using (135) in the same manner as above, we 
find 



*(s, a) = (I - o)-+ E (- iyk-i)(k-2)i2l 



/1oo^ (2m- l)\P,^(0)-h 

^^^^ + I s l^- cos e ' 

\h\<l, 0^a<l, -^<e<l, 

which differs from (84) only in the form of the remainder term. 

22. The integrals of Schaar and Landsberg. Let us first assume s real 
and positive and replace /(<, a) in (132) by the series (135), so that 

S \ 2^il2nTrir+~t) + 2niri{2n^^ ) ^~''^^- 
Since | 2mvi ±t\^ 2mr, the absolute value of the general term in the 
series does not exceed ( j^ — y + j^ — tj j e""', and the integral 
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being convergent, the same is the case with 

I V( ? I ! le-" di 

Jo ^ I V 2n^^^■(2n1r^ + 2ni:i{2mri — t) J \ ' 

By theorem XXII cor., we may therefore reverse the order of integration 
and summation in the expression for w(s, a), and write 

w /*« / g2niria g— 2tnria \ 

co(s, a) = - S J^ (^ 2niri(2nxi -TO "*" 2n7rt'(2nxi - j ^""'^'• 

Replacing i by 2nxi/s in each integral (this is a real substitution since s 
is real), we obtain 

since | s ± <i | S s, the absolute value of each integrand does not exceed 
(llnirs)e~^""', and the series 



n=I «/0 



1 °° 

rurs ~ 2x^s 1 n^ 



being convergent, the same is the case with the series 



;Jri Jo 1 2nir\s + ti^s -tij 



dt. 



By theorem XXII cor., we may therefore again reverse the order of 
summation and integration in w(s, a), which gives 

/»» 00 J /g— 2nir(t+o») g— 2nir(t— oi) \ 

and summing the infinite series 

(140) 

+ ^ilogj-;.-^(,^)d<. 

Passing to the case of s complex, we may now write for 9J(s) > 

'-(*' «) = 2^1 r+Ji ^°S 1 _ e-<'+^'^* 
(141) 

1 r* 1 1 

+ 2^ Jo s^-^TTi^^^l-e— <^'^'- 
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In fact, l/(s + ti) and l/(s — ti) are holomorphic in s for 9J(s) g «, 
I s I ^ A and i S 0; since | s ± ii | S e and 



log 



log 



J g— 2ir(«±oi) 

(which is seen by expansion in series, as in (106)), and furthermore 

Jr^i 1 

converges (the integrand becoming infinite as log tfort-^O and vanish- 
ing as e"^"' for i -» oo), both integrals in (141) converge uniformly for 
9i(s) ^ e, I s I ^ A. The appUcation of theorem XXIII now shows that 
the right side in (141) is a function of s holomorphic for 9I(s) > 0; the 
same being the case with w(s, a), and both sides in (141) being equal for s 
real and positive by (140), it follows that the equality subsists in the 
entire half-plane 9?(s) > 0. 

Writing 1 - e'^"'-''^'''^ = re**', we have r^ = 1 - e'^'* cos 2ira + e'*", 

sin 2wa 
6 = arctan 



and making 

Hi, a) = aT ^°S (^ ~ ^"^^ ^os 2ira + e'^"), 

(142) "" , . „ 

,, , 1 , sm 27ra 

x(f, o) = — J7- arctan -^z-, ^ — , 

^ 2ir e^"' — cos 2wa' 

we find by a simple algebraic transformation of (141) 

(143) „(..„). j;-^»)d,-jr-^).,, 9iw>o. 

Differentiating (141) in respect to x = 8t(s), we find 
(I - a) - - w*(s, a) = — ^ 

(144) = ~ 2^; X {s'+W ^""^ 1 - e-^''<'+<'*> '^^ 



2t Jo (s - <i)2 ^°S 1 - e-^'e--") '^'' 



the differentiation under the integral sign is legitimate by theorem XXI, 
since the integrals obtained converge uniformly for 3J(s) ^ e, \ s\ ^ A, 
as is seen in exactly the same way as for the corresponding integrals in 



120 T. H. GBONWALL. 

(141). Integrating by parts in the formula just obtained, we find for 
< o < 1 

~ 2^ J„ (s + uy ^^^ 1 _ e-2.(.+ao <^* = [2^- ■ M^ ^^^ l-e-^" ('+<"■) i=„ 

Jr>oo 2 g— 2)r(»+o») 

and the first term to the right equals 

e'"' — e~ 



l.log(l-e— 0=2^1og(2e''"'->'- 



oifat __ £> — vai ' 



27ris ^ ^ '^ 27ris ^ \ 2i J 

= 2:^ ["-(i - a)i + log (2 sin ira)]. 

The second integral in (144) is found by changing the sign of i, and we 
finally see that 

Jf»» 2 g— 2ir(t+ai) 

(145) 

/»» 1 g-2T((-oi) 

~'I g^r7^- l-e-'""-°^> ^^' ^W>0. 0<a<l, 
or writing 

, , cos 2x0 — e'^"' 

^^^^' ^^ ~ e^"' - 2 cos 2iro + e-^" ' 
(146) 

/, „N sm 2ira 

XiKh a) e2nt _2 cos 27ro + e-^"' ' 

and transforming (145) algebraically 

(147) u*{s,a)=j — 2 _|_ ^2 dt + j^ g2 _|_ ^2 dt, m{s) > 0. 

(147) is also seen to be true for a = by making a = in (143), differ- 
entiating in respect to x and integrating by parts. In the special cases 
a = and = 5 these expressions are considerably simplified, and we 
obtain 



«(«. 0) = - / ^^2 log i^z^2ir> dt, 
(148) I r s 1 



"*(s, I) = r 





2< d< 



S2 -|_ i2 g2,rt _^ 1 
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The integrals (148) are due to Schaar,*' and (143) and (147) to Landsberg.^* 
Using the identity 

1 '-^(- i)''-ifk-2 (_ l)m-l;2m-2 



s^ + f i£l s2* ^ s2'»-2(s2 + f2) 

in (143) and (147), and proceeding in the same manner as at the end of 
§ 15, there result asymptotic expansions which must coincide with (137) 
and (139). Comparing coefficients, we readily obtain 



(149) 



2*+l 

2 



'"dt. 



23. Kummer's trigonometric series. Let s be real and e ^ s ^ 1 — e; 
we then obtain from (128) 



, r(s 
log 



b=r['-^v^^+<-Hf- 



r(i 

and subtracting (125) multiplied by 2s — 1 

We now have the Fourier expansion 

z —^ h ^ , = -9 + Z^ (a„ cos 2TCirs + o„ sm 2nirs) 

uniformly convergent for ^ s ^ 1 and a fixed value of t, since the left 
side is a holomorphic function of s. It is seen from (34) that Oo = 
and that for n > 



*> Landsberg, G., Sur un nouveau d^veloppement de la fonction gamma, M6m. Ac. Belgique, 
vol. 55 (1897). 
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2 r g(2nT<— o« g-i+(2n7ri+()»"|»=i 2 f (2s l)e*"'^" "l'~' 

" 1 - e-' L2nirt - < ~ 2n7riT'r J,=„ "^ 1 + < L 2rm J.=o 

2 2 2i 



2n7rt — t 2mri + t nir(l + t) 
Sniri 2i 



so that 

(151) 

where 



4nV + (" W7r(l + ' 

sin 2nws 



re-.< _ e-(i-.)< 2s - 11 1 ^ ,,,sii 
L 1 - e- r- + T+T J • ? = L'^^^^ ~ 

^s ^ 1, < so 



nir ' 



_ 1 Z' 8"'^' 2 \ 2_ 

(152) c„W - ^ y^^2^2 _,_ ^2 1 + < j - 1 + < 



2t 



inV + f 



increases with n for any fixed positive value of t. We shall now prove that 
it is permissible to substitute (151) in (150) and integrate term by term, 
or in other words, that 



uniformly for e ^ s ^ 1 — e. Writing 



as n -^ 00 



„ -^ sin 2juirs 



M=i. ir/t 



so that 

sin 2j'ir3 



— Sy — S,+u 



PIT 

it follows from (38) that for e ^ s ^ 1 — e 

1 



Sy\^ 



vir sin ire 
and from the identity of summation by parts 

2^ Cy{Sy — Sy+i) = ^ (C„ — Cy^ijSy 4" C„On+l — Cn+p+lOft+p+l 
y=n+l i'=n+l 

we conclude, since Cy{t) — Cy+i{t) > for < > 0, that 
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(153) 



23 C.(0 -^ ^ ^ „:„ ^^ E (C.(0 - C,_i(0) • - 

MOl . li 

w + 1 n + p + 



I I Cn(0 I ■ I Cn+p+l(0 I "I 

"•"w+l "'"n + p + lj' 
From (152) we obtain 

(154) J%„mt = [log 4^^]^;" = 2 log 2n., 

and since Cn{t) < as i 5 4n27r2, 

I c„(0 \dt=\ Cnit)dt - c„(t)dt 

«/0 •-'0 t'4»2ir2 

= 2 log ( 1 + —r^) + 2 log 2n7r < 4 log 2nx, 
so that, by (153) 

Jo I y^+i VT T Sin ire L v=^+i J- \ v - IJ 

2 log 2n7r 2 log 2(n + p + l)ir "] 
+ n + l+ n + p + 1 J* 

Letting p -^ <» , and observing that log ( 1 + ~~zr\ ) *^ _ ■« and 

00 2 2 

23 ~7 ?\ = — , we finally see that 

1 r( V- /.N sin2j'irs'\ , I ^ 2 /I , 2log2nir\ . 

I I 2^ cXt) ]dt\< — : 1 - H n~i- ) -4 as n -» 00 

I Jo Vvin-i "T / I irsinireVn ' n + 1 / 

uniformly for e S s S 1 — e. Substituting (151) in (150) and integrating 
term by term with the aid of (154), we find for < s < 1 

, r(s) , ^,„ ,, ^ ^ log n + log 27r . - 

the series being uniformly convergent for e ^ s ^ 1 — e. Using (40) 
and (121) (the latter being identical with J 6), we now immediately 
obtain Kummer's series (113). 

Additional note. Application of complex integration to the theory of 
the Gamma function. We observed at the beginning of chapter III 
that when the Gamma function is defined by means of Euler's integral of 
the second kind, we are thereby restricted to the region of convergence 
of this integral, viz., the half plane 3J(s) > 0. This difficulty is overcome 
by the use of complex integration and other more advanced parts of the 
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theory of functions of a complex variable. The author's original inten- 
tion to devote a fourth chapter to an exposition of these methods had 
to be abandoned through lack of space, and it was found feasible only to 
give the following brief references to the literature. 

A paper by Birkhoff'^ contains what is probably the shortest way of 
arriving at the main properties of the Gamma function from the function 
theoretic point of view. Starting from the definition 

T{s) = lim V(l+ A + i)^ 



n-*« s{s + !)■ ■ -{s + n)' 
<p{s) = s*~i/^e~' V2ir, I arc s | < tt, 
which is essentially equivalent to (J 15), Birkhoff first proves the formula 
(J 14) log r(s) = (s - I) log s - s + log V2^ + u(s) 

with the Gudermann expansion of w(s) (J 13) and the fundamental 
property of this function, viz. | w(s) | < K/p, where X is a constant and 
p the distance of s from the nearest point on the negative real axis. Next, 
the formula r(s)r(l — s) = ir/sin tts (J 6), the infinite products for 
r(s) (J 3 and 5), and Euler's integral of the first kind are established, 
and after a digression on ^(s) = logs — w*(s), |w*(s)| < K/'p"^, Birk- 
hoff's paper closes with the expression of Euler's integral of the first kind 
in terms of the Gamma function (71 in the present paper). 

By the use of Cauchy's integral, a systematic study may be made of 
a whole class of summation formulas, containing that of Euler, and as an 
application, the formulas (76) and (82), Binet's integral (97) and those of 
Schaar (143, 147) may be derived from a common point of view, together 
with the various asymptotic expansions to which they lead. An excellent 
exposition of this aspect of the theory is given by Lindelof .'^ 

Attention is finally called to an important paper by Mellin,^^ in which 
is set forth the intimate connection of the Gamma function with a large 
class of other transcendental functions (for instance, the Riemann Zeta 
function). The main purpose of Mellin's paper is however to apply the 
Gamma function to the systematic solution of a class of linear difference 
equations (the simplest case of which is solved in / § 11) as well as to a 
corresponding class of linear differential equations, the coefficients of 
which are polynomials of the first degree. 

51 Birkhoff, G. D., Note on the Gamma function, Bulletin Am. Math. Soc, vol. 20 (1913-14), 
pp. 1-10. 

'^Lindelof, E., Le calcul des rfeidus, Paris, Gauthier-Villars, 1905 (Collection Borel). 

'' Mellin, H., Abriss einer einheitUchen Theorie der Gamma- und der hypergeometrischen 
Funktionen, Math. Annalen, vol. 68 (1910), pp. 305-337. 



